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ABSTRACT 

The Cluster Lensing And Supernova survey with Hubble (CLASH) is a 524-orbit multi-cycle treasury 
program to use the gravitational lensing properties of 25 galaxy clusters to accurately constrain their 
mass distributions. The survey, described in detail in this paper, will definitively establish the degree 
of concentration of dark matter in the cluster cores, a key prediction of structure formation models. 
The CLASH cluster sample is larger and less biased than current samples of space-based imaging 
studies of clusters to similar depth, as we have minimized lensing-based selection that favors systems 
with overly dense cores. Specifically, twenty CLASH clusters are solely X-ray selected. The X-ray 
selected clusters are massive {kT > 5 keV) and, in most cases, dynamically relaxed. Five additional 
clusters are included for their lensing strength {9Ein > 35" at = 2) to optimize the likelihood of 
finding highly magnified high-z {z > 7) galaxies. A total of 16 broadband filters, spanning the near-UV 
to near-IR, are employed for each 20-orbit campaign on each cluster. These data are used to measure 
precise (ct^ ^ 0.02(1 -fz)) photometric redshifts for newly discovered arcs. Observations of each cluster 
are spread over 8 epochs to enable a search for Type la supernovae at z > 1 to improve constraints on 
the time dependence of the dark energy equation of state and the evolution of supernovae. We present 
newly re-derived X-ray luminosities, temperatures, and Fe abundances for the CLASH clusters as well 
as a representative source list for MACS1149. 6-1-2223 (z = 0.544). 

Subject headings: cosmology: dark matter — cosmology: dark energy — gravitational lensing: strong 
— gravitational lensing: weak — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

The Universe has proven to be far more intriguing in 
its composition than we knew it to be even just 14 years 
ago. It is a "dark" Universe where ~ 23% of its mass- 
energy density is made up of weakly interacting (and, as 
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yet, undetected) non-baryonic particles (a.k.a. dark mat- 
ter) and ~ 73% is as yet unknown physics (a.k.a. dark 
energy) tha t is driving a n accele rated expansion of the 
metric (e.g., Komatsu et al. 2010; Riess et al.|2011 ). The 



Hubble Space Telescope [HST) has played a key role in 



both of these mysterious dark components (e.g.. 


Riess 


et al.||1998 


Perlmutter et al.||1999 


Clowe et al. |200e 


). 



high end of the cosmic mass power spectrum, provide a 
powerful way to constrain the frequency of high ampli- 
tude perturbations in the primordial density field. As 
such, they play a direct and fundamental role in testing 
cosmological models and in constraining the properties 
of dark matter (DM), providing unique and independent 
tests of any viable cosmology and structure formation 
scenario, and possible modifications of the laws of grav- 
ity. A key ingredient of such cluster-based cosmolog- 
ical tests is the mass distribution of clusters, both on 
(sub) Mpc scales and across the range of populations. 
The best and highest resolution maps of DM distribu- 
tion in massive galaxy clusters come from observations 
of strong gravitational lensi ng made by the Advanced 



Camera for Surveys (ACS, Ford et al. 2003) onboard 



HST. As part of a Guaranteed I'ime Observation (GTO) 
program, deep (20-orbit) multiband (4-6 filter) observa- 
tions were obtained for five galaxy clusters: Abell 1689 
( [Broadhurst et al.|[2005{ [Limousin et al.||2007[ |Coe efaL 



2 



Postman et al. 2011 



2010 



Zitrin et al. 



Abell 1703 ( Limo usin et al. 

2010^; Abell '2T 



OT7|; (JL002 ¥+TB54 ( |Jee et al |[2ro7| jZ^ ai.ii2009i 



20081 ISaha & Read) 
.'S^ (Eli'asdottir ct al. 



Umetsu et aL][2010| ); and MS 1358+6245 ( [Zitrin et al 
THese results have contributed to a reported ten- 
sion between observed and simulated galaxy cluster dark 
matter halos. Observed clusters appear to have denser 
cores than simulated clusters of similar total (virial) mass 



(e.g., Broa dhurst et al. 2008 Oguri et al. 2009 Sereno 



et al.|2UTU ). Understanding the true constraints from ob 



served concentration and mass profile measurements on 
A cold dark matter (ACDM) structure formation mod- 
els is one of the important problems that can be tackled 
with a deep, high angular resolution imaging survey of 
a significantly larger and more homogeneously selected 
sample of galaxy clusters. 

In May 2009, NASA successfully executed the final 
planned Hubble Servicing mission, SM4, w ith the instal- 
lation of the Wide Field Camera 3 ( WFC3 iKimble et aL 
2008 1 and the repair of A CS. Shortly thereafter, the Hub- 
ble Multi- Cycle Treasury (MCT) Program was conceived 
to permit ambitious programs (> 500 orbits) with broad 
scientific potential that could not be accomplished within 
the constraints of a single HST observing cycle and that 
would take full advantage of the final era of a newly re- 
furbished Hubble Space Telescope. 

The Cluster Lensing and Supernova survey with Hub- 
ble (CLASH) was one of three MCT programs selected. 
CLASH has four main science goals: 

1. Measure the profiles and substructures of dark mat- 
ter in galaxy clusters with unprecedented precision 
and resolution. 

2. Detect Type la supcrnovac out to redshift z ^ 2.5 
to measure the time dependence of the dark energy 
equation of state and potential evolutionary effects 
in the supernovae themselves. 

3. Detect and characterize some of the most distant 
galaxies yet discovered at z > 7. 

4. Study the internal structure and evolution of the 
galaxies in and behind these clusters. 

To accomplish these objectives, the CLASH program 
targets 25 massive galaxy clusters and will image each 
in 16 passbands using WFC3/UVIS, WFC3/IR and 
ACS/WFC. CLASH has been allocated 524 T/ST orbits, 
spread out over Cycles 18, 19, and 20. The majority of 
these orbits (474) are for cluster imaging and, simulta- 
neously, for the parallel SN search program. An addi- 
tional 50 orbits were allocated as a reserve for supernova 
follow-up observations. Based on a current census of the 
HST data archive, CLASH will produce a six-fold in- 
crease the number of lensing clusters observed to a depth 
of 20 orbits and, more importantly, will vastly increase 
the number lensing clusters with extensive multi-band 
HST imaging. 

Motivations for each of th e ma in CL ASH science goals 
are provided in §2.1 [ \ 2.2[ j |2.3[ and ^.4| Subsequent 
sections describe the cluster sample (§31), survey design 
(Q, data pipeline (^|, and supporting observations us- 
ing other facilities (§OT. Data products intended for pub- 
lic distribution to the community are briefly described in 



S|7] AB magnitudes are used throughout (Oke 19741. The 
cosmological parameters H^ — lOO/i km s~^ Mpc~^, h = 
0.7, = 0.30, A — 0.70 are assumed in this paper. 

2. SCIENTIFIC MOTIVATION 

2.1. Galaxy Cluster Dark Matter Profiles and 
Formation Times 

Recent observations suggest that galaxy clusters 
formed earlier in our universe than in simulated ACDM 
universes. These observations include the detection of 
perhaps unexpectedly massive galaxy clusters at z > 1 



(Stanford et al. 20061 Eisenhardt et al. 


2008 


Jee et al. 


20091 IHuang et al. 20091 Kosati et al. 




IPapovich 


et al.pOlO^ 'Schwope et al.''2010| 


Gobat et al. |2011 Jee 


et al.|2011; Foley ct al. 2011) and 


the finding that some 



cores t han clusters of similar mass prod uced in simu- 
lations dBroadhurst ct al. 2008 



Broadhu rst &: Barkana 
20081 lUguri et al.|2009iiKichard et al.,20lOl ISereno et al. 



2010| |Zitrin et al.||20Ha[ ). While the evidence to date 
tor early cluster growth is suggestive, possible explana- 
tions include departures from Gaussian initial density 
fluctuation spectrum or higher levels of dark energy in 



the past, so-calle d Early Dar k Energy (E DE) (Fedeh fc: 
Bartclniann 2 0071 ISadeh fc R cphacli 2008] [Fraricis et aT 
2009[ jCrossi fc Springe! 2009j). If a significant quan- 
titv 01 EDE (tor example S2n; 



tity ol EDE (tor example Ude 0.1 at z = 6) sup- 
pressed structure growth in the early universe, then clus- 
ters would have had to start forming sooner to yield the 
numbers we observe today. These scenarios remain al- 
lowable within current observational constraints as de- 
scribed in the above papers, although some non-Gaussian 
models can be ruled ou t by using the cosmic X-ray back- 
ground measurements ( Lemze et al. 112009) . The CLASH 
data permits significant advances to be made towards 
supporting or rejecting observational evidence for early 
cluster growth by measuring core densities for a larger, 
less biased sample of clusters. 

In cosmological simulations, CDM-dominated halos of 
all masses consistently evolve to have a roughly "univer- 
sal" density profile that steepens with radius. Functional 
for ms that fit such a profile wel l include the "NFW" pro- 
file jNavaiTO~erarn^^ and the Einasto / Sersic 



pro! 



2010 



IS re 



Scrsic 1963; Einasto 19651 iNavarro et al.ll2004 
Furthermore, each simulated halo's core density 
ated to the background density of the universe at 
the halo's formation time. Halos that form later, includ- 
ing the most massive galaxy clusters, are found to have 
the least dense cores in a relative sense. Determining 
the relationship between the shape and depth of a halo's 
gravitational potential and its total mass as a function of 
time thus provides fundamental constraints on structure 
formation. 

In practice, the relative core densities, or "concentra- 
tions", are measured (both in simulated and observed 
halos) as Cyir = ryir/r^2, a ratio between the virial ra- 
dius and the inner radius at which the density slope of the 
fitted profile is isothermal (p oc r^^). Analyses of grav- 
itational lensing spanning such a large range of radius 
allows one to map the (primarily dark) matter profiles of 
observed halos and measure their concentrations. 

Dark matter profiles are best mapped in cluster cores 
using strong lensing analysis of multiband HST imaging. 
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Table 1 

Concentration Measurements for Previously Well-Studied Lensing-selected Clusters 

















Constraints^ 


Publication 


Cluster 


z 


[lO^^MQh-'^] 




xVd.o.f. 



SL-l-WL-l-mag 

SL-l-WL-l-mag 

SL-l-WL-l-mag 

SL-l-WL-l-mag 

SL-l-WL-l-mag 

RE-I-WL 

RE-I-WL 

RE-I-WL 



Umetsu et al. 



Umetsu et al. 



1 2011b i'' 

|2011b|fa 

Umetsu e t al. ( 20lTb i'' 
20lTb i'' 



Umetsu et al. 



Umetsu et al. 



Oguri et al. 



Oguri et al. 



Oguri et al. 



I 2009 1 



2011b !*> 



I 2009 1 



I 2009 1 



A1689 0.187 1.34 

A1703 0.281 1.29 

A370 0.375 2.26 

CL0024-I-17 0.395 1.37 

RXJ1347.5-1145 0.451 1.73 

SDSS J1531-I-3414 0.335 0.7; 

SDSS J1446-I-3032 0.464 0.8 

SDSS J2111-0115 0.637 0.9 



+0.20 
0.16 
+0.22 
0.19 
+0.26 
0.23 
+0.20 
0.18 
+0.12 
0.11 
+0.29 
0.24 
+0.3 

0.22 
+0.41 
0.32 



,+ 1.72 
-1.62 
3+1.04 
="-0.91 

4.56 ±0.33 

rj y„ + 0.97 

'■' '-0.87 
"+0.37 
-0.35 
3+3. 
-1.5 
,+3.9 



13.82: 
6.1 



5.96! 
7.9j 
8.3"* 



4.73/17 
7.14/19 
14.07/24 
11.47/20 
45.06/25 
8.1/6 
6.4/6 
7.5/6 



Note. — Of these 8 clusters, only one - RXJ1347. 5-1145 - is in the CLASH sample. 
^ SL — strong lensing; WL — weak lensing; mag — magnification b ias (number counts); R E — Einstein radius 
^ The parameters here arc from standard NFW fits to profiles in |Umetsu et al.| ( |2011b[ l, and not from the generalized 
NEW (gNFW) fits given in that paper. 



©' 




1=1 



(C 10' 



o 

US 
U 

3 10' 



CO 



10 



A1689 (z=0.183) 
A1703 (2=0.281) 
A370 (2=0.375) 
C10024+1654 (z=0.395) 
NFW fit (stacked 4 clusters) 
... I i_ 

1 00 1 000 

R [kpc//i] 



6 

'u 



■ Hubble strong lensing 

■ Subaru weak lensing 
- Combined 




Figure 1. Mass profiles measured for four well-studied, strongly 
lensing galaxy clusters that are not included in the CLASH sam- 
ple. All have similar mass profiles as measured from Hubble obser- 
vations of strong- lensing and Subaru o bservations of v^eak- lensing 
distortion and magnification | |Ume tsul;t al. 2011b. their Figure 6). 
The averaged mass profil e is in r emarkably good agreement with 
the standard NFW form ([Umetsu et al. 2011a( their Figure 1) as 
shown by the gray area (2-cr confidence interval of the NFW fit), 
though with a higher concentration than predicted from cosmolog- 
ical simulations. Both strong and weak lensing probes are required 
to map the continuously steepening mass profile from the inner 
core (~ 10 kpc//i) out to beyond the virial radius (~ 2 Mpc/h). 

The lensing-based mass profile mapping is extended to 
the virial radii of clusters using weak lensing analysis 
of wider field ground-based multiband imaging, suc h as 
from Subaru. Results from Umetsu et al. ( 2011b|a ) for 
four of the currently best-studied (non-CLASH) clusters 
are shown in Figure [l] 

Joint modeling of strong and weak lensing (SL-I-WL) 
yields significantly better constraints on co ncentrations 
than either probe alone. Quantitatively, |Meneghetti| 



5 10 15 

Concentration, c-=r-/r^ 

* vir vir' s 



20 



Figure 2. Joint strong and weak lensing analyses are required 
to obtain tight constraints on cluster concentrations as shown here 
for CL0024+17, a non-CLASH cluster (Umetsu et al. 2010, from 
their Figures 15 and 17). Confidence levels of 68.3%, 95.4%, and 
99.7% are plotted in the Cvir^Afyir plane. 

have been found to have overly high concentrations 
(dense cores) compared to haloes in N-body simulations 
with similar masses, as shown in Figure ( Broadhurst 
et al.|2008| [Oguri et al.|2009[|Sereno e t al. 2010) andde- 
tailed in 'i'able[lj Recent siniulations (Prada et al.j2011 1 
yield cluster concentrations th at are over 50% higher 
than previous simulations (e.g., [Duffy et al.||2008[ ). This 
is a product of upturns as a function of both mass and 
redshift found in these newer simulations, which are not 
yet understood. Similarly, clusters have also been found 
to have somewhat larger than expected Einstein radii, 
a direct and particularly a ccurate measure of the pro 



jected mass in a halo's core (Broadhurst &: Barkana|2008 
Richard et c[r][2010[ [ZiFrin et al.[[20fTa|E| 



et al. (2010a) found their joint SL-j-WL analyses of sim- 
ulated clusters yield concentrations to ^ 11% accuracy, 
while WL-only and SL-only analysesyielded ^ 33% and 
~ 59% scatters, respectively. Figure [2] demonstrates how 
SL and WL analyses combine to yield robust constraints 

on the mass and concentration of CL0024-fl7 (Umetsu the fine of sight (iHennawi e t al. 2007 



et al. 2010). 



Unfortunately, the best studied clusters to date have 
also been among the strongest gravitational lenses 
known. Such lensing-selected clusters are highly biased 
toward halos with high concentrations, both intrinsically 
and as projected on the sky d ue to halo elongation along 



ford 2009 Meneg: 



letti et al 



i'he best studied (SL-fWL) galaxy clusters to date are estimated to lead to 



Oguri & Bland 
2010bn2Un| ). These biases' 
systematically higher concen- 
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Postman et al. 2011 
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Figure 3. Mass profiles of the best studied clusters to date are 
revealed to have higher central density concentrations than sim- 
ulated clusters of similar mass and redshift. Reconciliation may 
be within reach given results from the latest simulations and an 
estimated lensing bias which CLASH will avoid. The plotted lines 
are mean concentrations at three different redshifts for clusters in 
these simulations as calculated from the fitting formulae provided 
in those papers. However note that halos of this great mass are 
rare (or even non-existent at these redshifts) in these simulations, 
so these results are mainly extrapolations, as designated by the 
dashed lines. The thinner da shed lines above ill ustrate a 50% ob- 
servational bias applied to the |Prada et al.|||2011[ l results. This bias 
has been roughly estimated for non-UijABrl clusters such as these 
which were selected for study based on exceptional lensing s trength 
^ennawi et al. 2007, Oguri & Blandford 2009, Meneghet ti et al.| 
l^iUa 2Uii). Kesults from the currently best-studied clusters are 
Wotted here as squares (Umetsu et al. 2011b) and circles (Oguri 
[et al.|2009i) labeled with abbreviated names and described further 



alj2009|) 
'i'at)le[Tl 



trations by as much as 50% or more. However, such a 
bias is insufRcient to account for the discrepancy between 
the observations and the predictions that is until a very 
recent analysis of new simulations was performed (see 
Figure [s] and discussion below) . 

More robust conclusions require analysis of a larger, 
unbiased cluster sample. Progress toward this goal has 
been m ade by LoCuSS, th e Local Cluster Substructure 
Survey (Smith et al. 20051. A large sample of 165 clus- 
ters between 0.15 < z < 0.30 was selected based on X- 
ray brightness. Strong lensing analyses of 20 of these 
based on HST imaging (mostly single -ban d "snapshots' 



were pr esented by Richard et al. (2010[ ). Okabe et al. 
(2010a I published weak lensing analyses ot 30 LoUuyS 
clusters, 22 of these being more "secure" based on multi- 
ban d Subaru imaging, an d 9 of these 22 overlapping with 
the Richard et al.| (2010) subset. Subaru images are es- 
peciaily desirable for WL studies as they enable excellent 
galaxy shape measurements to be performed over a wide 
area. Multiband imaging is also critical to properly select 
background galaxies and avoid significantly diluting the 
weak lensing signal (and thus the virial mass and con 



centration) with unle nsed foreground galaxies ( Medezin- 
ski et al. 112007 2010). Stacked WL-only analys es have 



been performed on large cluste r samples ( (Johnston et al. 



2007 



Mandelbaum et al.|2008 1 , but we re-emphasize the 



need for combining strong -|- weak lensing analyses in ad- 
dition to cross-comparisons with mass profile estimates 
from other techniques. For example, mass concentra- 



Ettori et al. 2010), although these are subject to 



uncertainties due to assumptions about hydrostatic equi- 
librium. Importantly, cluster elongation along the line of 
sight (a potential bias in concentration measurements) 
can be meas ured by the com bination of l ensing and X - 



ray analysis (Morandi et al. 



The caustic te chnique (IDiaterio 



2010 



3. N ewman et al.||2011 1 
¥^ellcr 1997riRines L 



Diaferio 20061 is, like fensing-based methods, indepen- 
dent of the dynamical state of the cluster and also pro- 
vides an important cross check on mass estimates from 
lensing and gas kinematics. Further discussion of some 
of the previous results from var i ous method s are given in 
iComerford fc Nataraja nl (|2007[ ); [Coe| ( |2010[ ); fRines et al ' 
l^OlO); King & Mead (?OTTf~^ 

Reconciliation ol high observed concentrations with re- 
sults from simulations may ultimately come from sig- 
nificantly reducing the observational sample bias along 
with finding higher concentrations in simulated clusters. 
Baryons, for example, are currently absent from those 
cosmological simulations large enough to produce mas- 
sive clusters. Baryons can result in significant "adiabatic 
contraction" on galaxy scales, however they constitute a 
much smaller fraction of the mass on cluster scales. Sim- 
ulations including baryons show that cluster halo con- 
centrations are likely only varied by 10% or so relative 
to DM-only halos, and that the direction of variation (in- 
crease or decrease) i s not even clear, depending on the 
gas physics assumed ( [Duffy et al.| 2010||Mead et al.|20l"o{ 
King & Mead 20 11| ). Nonetheless, measuring tfKTveloc 



ity dispersion of the Brightest Cluster Galaxy (BCG) as 
an additional constraint on the inner (< 80 kpc) mass 
profile, where its stellar mass is a non- negligible com- 
ponent of the matter distribution, provides for a more 
thorough mapping of the total mass profile. Some clus- 
ters have been shown t o have inner mass profiles that are 



shallower than NF W (Sand et al. 2004 2008 Newman 



et al. 2009 2011). Deviations of cluster mass profiles 



from iNl^'W or Einasto forms at small and/or l arge radii 



may slightly b ias concentration measurements (Oguri & 
Hamana|[2011[) . The use of multiple probes of the mat 



ter distribution, as CLASH is designed to do, will enable 
such biases to be measured and the significance of devi- 
ations determined. 



tions can be measured from X-ray profiles (Buote et al 



More recently, an analysis by Prada et al. (2011 1 found 
that clusters in the Bolshoi and MultiDark simulations 
have concentrations ^ 50% higher than clusters in previ- 
ous simulations. While the robustness of this new result 
is still being assessed, it raises the possibility that the 
combination of new observations of an unbiased sample 
of clusters and new simulations may be able to bridge 
the concentration gap. We stress here that estimates 
of the observational bias from previous cluster studies 
have large uncertainties and likely varies for each cluster. 
CLASH will determine mass profiles and concentrations 
for a new cluster sample free of lensing selection bias. 
As we demonstrate in p.2[ CLASH is designed to detect 
(or rule out) with 99% statistical confidence average de- 
viations of 15% or more from predicted concentrations. 
However, given the outstanding uncertainties in the ex- 
pected concentrations, we prefer to recast the problem 
as follows: CLASH will deliver robust observational con- 
centration measurements for a sample of clusters that 
simulations will be tasked to reproduce. Ultimately this 
will lead to a better calibration of many mass estimation 
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techniques and, consequently, to a better understanding 
of structure formation on cluster scales and perhaps of 
our cosmological model as well. 

2.2. Improved Constraints on the Dark Energy 
Equation of State and SNe Evolution 

The biggest cosmological surprise in decades came from 
observations of high-redshift type-la supernovae (SNe 
la), providing the first evidence that the e xpansion of 
the Universe now appears to b e accelerating (iRiess et al. 
1998[ Perlmutter et al. 19991, and indicating the Uni- 
verse is dominated by darJt energy." The presence of 
dark energy has galvanized cosmologists as they seek to 
understand it. Observations of high-redshift SNe la have 
continued to lead the way in measur ing the properties of 
dark energy (e.g., Riess et al. 2011). The goal for cos- 
mologists now is to measure the equation of state of dark 
energy, w = P/{pc^), and its time variation in the hope 
of discriminating between viable explanations. A depar- 
ture of the present equation of state, wq, from —1 or 
a detection of its variation, dw/dz^ would invalidate an 
innate vacuum energy (i.e., the cosmological constant) 
as the source of dark energy and would point towards a 
present epoch of "weak inflation." A difference between 
the expansion history and the growth history of structure 
expected for w{z) would point towards a breakdown in 
General Relativity as the cosmic scale factor approaches 
unity. 

HST paired with ACS is a unique tool in this inves- 
tigation, providing the only means to collect SNe la at 
1 < z < 1.5, which, in turn, provide the only constraints 
we have to date on the time variatio n of w. From the 
23 SN e la at z > 1 with HST data ( [Riess et al^|2004 
20071 we have learned: 1) that cosmic expansion was 



once decelerating before it recently began accelerating 

2) that dark energy, i.e., an energy density with w < Q, 
was already present during this prior decelerating phase, 

3) that SNe la at a look-back time of 10 Gyr appear both 
spectroscopically and photometrically similar to those 
seen locally and 4) no rapid change is seen in w{z) and 
thus no departure is yet seen from the cosmological con- 
stant, though the constraint on the time variation re- 
mains an order of magnitude worse than on the wq. 

SNe la play a central role not only as distance indica- 
tors for cosmography, but also as major contributors to 
cosmic metal production and distribution. The measure- 
ment of high-redshift SN la rates is therefore integral to 
understanding the history of chemical enrichment. The 
high-redshift rate cannot easily be predicted from the 
star formation rate because the nature and, hence, the 
timescales of the process behind the growth of the white 
dwarf toward the Chandrasekhar mass is not known. The 
two leading competing scenarios are accretion from a 
close binary companion - the single- degenerate scenario 
( [Whelan fc Iben|1973t |Nomoto| 19821) or merger with an- 
other white dwarf, following loss of orbital energy and 
angular momentum by emission of gravitational waves 



the do uble-degenerate scenario fiben & Tutukovl 



1984 



Webbink^ 1984) . One way to constrain the ditterent pro- 
genitor scenarios is to measure the delay-time distribu- 
tion (DTD) of SNe la. This is the distribution of times 
that elapse between a brief burst of star formation, and 
the subsequent SN la explosions. Observations h ave sug- 
gested various different forms for the DTD (e.g., (Dahlen 
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Figure 4. Dark Energy and Evolution Sensitivity. Three models 
consistent with current data are shown: wq = 0.9, = (black); 
uiQ = 1, Wa = 0.8 (red); and = l,Wa = (dashed). Also plotted 
blue is the SNe la evolution model of iDominguez et al.| | |2001| , 
where the peak luminosity changes 3% per solar mass cliange m 
the donor star. The error bars show the constraints at present, as 
projected after CLASH, and if HST continues to collect SNe la at 
the present rate for 7 years. 




eTaTl [20041 |2008l [Mannucci et~ar[[2006l [Pritchet et all 
2008[ )). However, a number of more recent measurements 
and analys es point to a DTD that is a pow-er law of in- 
dex w -1 ([Totani et al.[[2008[ [Maoz et al.[[20Tol [Maoz 



fc Badenes'2010; "-Brand t et al.i20 10','Hori uchiTljeacom 
2010 Maoz et al. 2011; Graur et al. 201 1|. Specihcally, 
Graur et al. (2011|) have shown that such a DTD fits well 
the measured SiN la rate out to z « 2. Small sample sizes 
are the major limiting factor at high redshifts. Large 
high-z SN samples are therefore needed to resolve the 
issue, and control possible biases in cosmological studies 
(due to the evolving SN la channel mix). 

The CLASH survey uses ACS in parallel with the 
cluster program to continue the discovery of SNe la at 
1 < z < 1.5, the objects which tell us about the variation 
in w. With WFC3 in parallel, CLASH w ill yield SNe la 



at 1.5 < z < 2.5 (Rodney et al. 2011). Observations 



in this fully matter dominated epoch provide the unique 
chance to test SN la distance measurements for the dele- 
terious effects of evolution indepen dent of our ignorance 
of dark energy ( Riess fc Livio|2006 ). Because the SNe la 
are detected when these cameras are in parallel, they are 
far from the cluster core (^ 2 Mpc at the median cluster 
redshift of z = 0.4) and, hence, the effects of lensing are 
small (and correctable), making the SNe usable for im- 
proving the limits on the redshift variation of the dark 
energy equation of state. At z < 1, SN la distance mea- 
surements are most sensitive to the static component of 
dark energy, wq. At 1 < z < 1.5, the measurements 
are most sensitive to the dynamic component, Wa- By 
z > 1.5, the measurements are most sensitive to evolu- 
tion if present (e.g., the changing C/0 ratio of the donor 
star), providing the means to diagnose and calibrate the 
degree of SN la evolution in dark energy measurements. 
Figure [4| shows how variations in the DE equation of state 
or an evolution in the white dwarf C/0 ratio can change 
the SN la distance modulus as a function of redshift. 
Accounting for the systematic unc ertainty introduced 



by the cosmic star formation history, Graur et al. (2011 ) 
predicted the SN la rate out to higher redshilts (their 
figure 13), which we use here, in conjunction with the 
CLASH observational parameters, to estimate the num- 
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Table 2 

Estimated Number of z > 1 SNe la to 
be Discovered During the CLASH 
Program 



Redshift 


Estimated Number 


Range 


of SNe la Found 


1.0< z < 1.5 


7-11 


1.5< z < 2.0 


4-13 


2.0 < 2 < 2.7 


0-4 


1.0 < 2 < 2.7 


11 - 28 







Note. — These estimates arc the 68% 
confidenee intervals, aceounting for both 
statistical and systematic uncert ainties. 
The assume d SN la rate is from IGraurl 
|et al.|j2011t . 

ber of z > 1 SNe la that will be discovered over the 
course of this program. These estimates are presented in 
Table [2] and are the 68% confidence intervals, accounting 
for both statistical and systematic uncertainties. 

2.3. Detection and Characterization of z > 7 Galaxies 

One of the most important goals in observational cos- 
mology is to find the first generation of galaxies and un- 
derstand their roles in the reionization of the intergalac- 
tic medium (IGM) between z ~ 6 to z ~ 12. Substantial 
progress has been made in recent years in finding ob- 
jects at z 6-7 (Bouwcns et al. 2006, 2009; Bunker' 
et al.|20101|Oesc h et ai. 2010a b, McLure et al. 2010, Ca-i 



pak et al |[201ip , which" are, at'most, ~ 700 - 900 Myi 
oi age. We wish to find objects at even earlier epochs, 
in order to understand (1) how the first galaxies were 
formed through a hierarchical merging process; (2) how 
the chemical elements were generated and redistributed 
through the galaxies; (3) how the central black holes ex- 
erted influence over the galaxy formation; and (4) how 
these objects contributed to the end of the "dark ages." 

The majority of galaxies at z = 6 — 7 have been dis- 
covered via two methods: (1) deep pencil-beam s urveys, 
such as th e Hub ble Uhra Deep Field (HUDF; Beck-' 
with et al. (2006)) and t he Great Observa t ories Origins 
Deep Surwy(ITOODS; [Giayalisco eFaT] ( [20041 )); and 
(2) degree-size surveys with lOm-class ground-based tele- 
scopes, such as the Subaru Deep Field. Both use a color 
selection to search for "dropout" candidates - galaxies 
with deep IGM absorption at the wavelengths shortward 
of the redshifted Lya break. These efforts have proven 
successful at finding objects at z < 7, but there is grow- 
ing evidence for a rapid d ecrease in the number of candi 
dates at higher redshifts (|Iye et ah 2006 Bouwens et al 



2006 
73T 



Castellano et al 
7 and 59 z 
These results are based on 



2010| . H S"1 V WFCS data yieH s 
8 candidates ( Bouwens et al.|[2010 ). 
led on ~ 500 orbits of MUDl" obser- 
vations. In comparison, > 600 galaxies are found at z ^ 6 
(iBouwens et al. 1 120071 iSu et aL||2011|). To date, there 



have been only tw o or three spectroscopically confirmed 
galaxies at z > 7 ([Lehnert et al. 2010 Ono et al. 2012 



Schenker et al. 2012 ) a nd one S]3ectrum tor a Gamma 
Kay Burst at z ^ 8.2 (Salvaterra et al. 2009 Tanvii 
et al. 2009). The reason, in part, is a lack ot photons 
at redshift z = 8 and 10, an L* galaxy would have an 
apparent magnitude in the first NIR detection band of 
28.2 and 29.6, respectively. 




25 25 27 28 

Magnitude (AB, detection band) 



29 



Figure 5. The estimated number counts in the total CLASH 
survey at 2 ^ 7.5 and 10 are plotted as the solid "lensed" blue and 
red lines, respectively, as functions of magnitudes in the detection 
bands FllOW and F160W. Note the large uncertainties (shaded re- 
gions). A NIR survey of a comparable area (100 square arcmin) in 
unlensed fields would hav e a ~ lOx lower effi ciency at these magni- 
tudes (in agreement withlMaizy et al. 12010!). We have assumed an 
evolving Schechter lumi nosity function with dM* /dz = 0.36 as de- 
rived by (jBouwens et al. 2006 2008 ) and plotted as the da shed lines 
with shaded uncertainties. Observed number counts from lBouwensI 
et al. ( 2009 ) are plotted as squares. Lensing estimates were derived 
using the Zi trin et al.| ( |2009i ) CL0024+17 mass model to simulate 
the magnifications ot sources and reduction of source area. These 
estimates are in agreem ent with observed len sed counts in cluster 
fields plotted as circles jBouwens et al.|2009| . 

Gravitational lensing by clusters amplifies the fiux of 
background sources considerably. These cosmic tele- 
scopes improve the efficienc y of searches for relatively 
bright h igh redshift galaxies ( Bouwens et al.|2d09 Maizy 
et al.|2010| ). In fact, the majority ot m < 25.5 AB, z > 6.5 



galaxy candidates have been found in cluster fields (e.g 



kneib et al.|2004[ [Br adley e t al.||20q8l [Zheng et al.||2009, 
Bradley et al.||2011[ 'Schenk eret'al.||2012p . Interestingly 
most ot these candidates are found in regions with ampli- 
fication of ^ 10. Furthermore, lens models can be used to 
help discriminate between highly-reddened objects and 
truly distant, high redshift objects, as the projected po- 
sitions of the lensed images are strong functions of the 
source redshifts. 

Luminous z > 7 galaxies are extremely valuable as 
their spectra can be used to determine the epoch of the 
IGM reionization. This is because only a tiny fraction of 
neutral hydrogen is needed to produce the high opacity of 
Lya observed at z ~ 6. The damped Lya absor ption pro- 
file that resul ts from a completely neutral IGM ( Miralda- 
Escude||1998l ) can be measured even at low-spectral res- 
olution. I'urthermore, direct measurements of the early 



star-f ormation rate (via Lya and Ha emission, lye et al. 
2006 1 can be derived from the spectra of bright high-z 
galaxies. CLASH may detect dozens of relatively bright 
(magnified to m < 26.7 AB) z > 7 galaxies, including 
some bright enough for spectroscopic follow-up. A NIR 
survey of a comparable area (100 square arcmin) in un- 
lensed fields would have a '-^ lOx lower efficiency. An 
estimate, albeit a highly uncertain one, of the CLASH 
high-z detection efficiency enhancement is shown in Fig- 
ure [5l 

2.4. Galaxy Evolution 

ACDM provides a robust theoretical framework for the 
evolution of dark matter haloes. However, the forma- 
tion of galaxies within these haloes is governed by com- 
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plicated baryonic interactions that are difficult to simu- 
late on cosmological scales. The gap in our understand- 
ing of the connection between the formation of galaxies 
and that of their parent halo es is illustrated by observa- 



tions of galaxy "downsizing" (Cowie et al. 11996). ACDM 



predicts that structures form hierarchically, with small 
haloes forming early and later assembling into larger 
haloes. In contrast, massive galaxies appear to have 
formed at early times, as evidenced by the redshift at 
which their star formation peaks (e.g. 
Guzman et al.| 1997[ |Brinchmann 



Cowic ct al."1996t 
Mis 2U00; J uneau 
et al.||2005[ ) an d the old ages of their stellar populations 



at z 



UTe.g.; iFaber et al.||1995[ |Worthey ][T996, ,Proc- 

i'he stellar 



tor fc Sansomj 2002 ; T homas et ai.||2005p 



mass density of the universe is also increasing with cos- 
mic time, as galaxies form stars, then shut down their 
star forn iation and, ultin iately, accum ulate on the red 
sequence pell et al.[2004[ |Faber et al.||2007[ ) . 

Several open questions remain about the origin and 
nature of this stellar mass growth. Is it dominated by in 
situ star formation, or through mergers of existing stellar 
systems? How does the balance of these two processes 
change with galaxy mass and cosmic time? In a system 
dominated by hierarchical assembly, this is akin to asking 
whether most of the mass acc reted by gal axies comes in 
as gas (e.g., the "cold flows" of Keres et al. 2005 ; Dckcl & 
Birnbo im'2006!) o r as stars (e.g., "d ry mergers" , |Khochf"ar 
&: Hurk crt 2003) [Bell et al ||2004[ ). Theoretical expecta- 



tions predict that the in situ population should be very 
centrally concentrated, while an accre ted stellar compo - 
nent should extend to larger radius (Oser et al. 2010). 
Thus massive galaxies are predicted to have substantial 
gradients in the origin of their stars, with the innermost 
stars having formed in situ and the outer stars largely 
accreted through merging. 

Recent work has demonstrated t hat massive, passive 
galaxies already e xist at z ^ 2 (e.g., Trujillo et al.||2006 
Kriek et al.|2009 ) and that these systems are much more 



compact than local ma s sive galaxies (Trujill o et al.|2006 



|Longhetti et al.||2007| |van Dokkumj |2008j) , suggesting 
that they must grow significantly in size Between z ~ 2 
and z ^ 0. This process likely happens through "dry" 
(dissipationless) merging, since such massive galaxies are 
observed to have old stellar populations locally. There 
are also substantial populations of massive star-forming 
galaxies over a similar redshift range. T hese galaxies 



tend to exhibit large star forming "cl umps" (ICowie et al. 
1995| |van deiTB crgh 1996; Elmcgrce n et al.|2DIMa|bp , al- 



though at least some o f these galaxies m ay torrn disks 
with coherent rotation ( Genzel et al.|[20T0 ). 

There are several observational challenges to measur- 
ing and characterizing these two modes of growth. It 
is difficult to resolve substructure in high-redshift galax- 
ies, even from space. Both resolution and signal-to-noise 
considerations limit observations to the upper end of 
the mass function. Searches that use a small number 
of imaging filters struggle to identify the redshifts (and 
therefore the masses) of the targets and can only poorly 
constrain their spectral energy distributions (SEDs) and 
stellar populations. Studies that rely on emission lines 
probe only the assembly of star-forming galaxies. 

The CLASH program is well suited for studying galaxy 
assembly at 2; < 3. The high magnification and spa- 
tial stretching of strongly lensed distant galaxies, coupled 



with the broad 16-band photometry and resulting pho- 
tometric redshifts, enables the measurement of star for- 
mation rates and stellar ages in each lensed galaxy over 
many resolution elements. The magnification makes it 
possible to resolve substructures in high-redshift lensed 
galaxies further down the galaxy mass function than is 
possible in unlensed galaxies at similar redshifts. In ad- 
dition, with CLASH, we can contrast the relative prop- 
erties of galaxy cores versus their outer regions, and 
(where they exist) the properties of substructures such 
as clumps, spiral arms, and ongoing mergers. In the full 
CLASH sample of 25 clusters, we expect to find several 
strongly lensed background galaxies per cluster that are 
suitable for such analyses, providing a sample of 50 to 
100 galaxies with 1 < z < 3. 

3. CLASH CLUSTER SAMPLE 

The CLASH program is robustly measuring galaxy 
cluster dark matter profiles and concentrations for a sys- 
tematic comparison with those realized in cosmological 
simulations. Specifically, our cluster sample size and se- 
lection criteria were chosen to allow the robust measure- 
ment of deviations from the predicted cluster concentra- 
tion distribution of ^ 15% or more at high statistical 
confidence (~ 99% C.L. ) given a relatively unbiased en- 
semble of clusters {{2.1). 



3.1. Cluster Sample Selection 

To date, robust joint SL-I-WL analyses have only been 
performed on a small, highly-biased sample of five to ten 
clusters (Table [T] and Figure [s]). These clusters were se- 
lected for study primarily based on their gravitational 
lensing strength, which tends to preferentially select sys- 
tems with higher central matter concentrations. To es- 
tablish a sample that is largely free of lensing bias, we 
selected 20 massive clusters from X-ray-based compila- 
tions of dynamically relaxed systems. Sixteen o f these 20 
clusters were taken from the Allen et al. ( 2008 ) compila- 
tion of massive relaxed clusters. Clusters m our X-ray se- 
lected subsample all have Tx > 5 keV and exhibit a high 
degree of dynamical relaxation as evidenced by Chandra 
X-ray Observatory images that show well-defined central 
surface brightness peaks and nearly concentric isophotes. 
The clusters, in general, also show minimal evidence for 
departures from hydrostatic equilibrium in X-ray pres- 
sure maps. Most of the 20 clusters have smooth and 
only mildly elliptical ((e) = 0.19) X-ray emission and a 
BCG within a projected distance of 23 kpc of the X-ray 
centroid. No lensing information was used a priori to 
select them in order to ensure that they are not prefer- 
entially aligned along the line of sight, in contrast with 
a purely lensing-selected sample, where the surface mass 
density is, on average, biase d upwards along the hue of 
sight by intrinsic triaxiality ( Hennawi et al. 2007 Oguri 
et al.||2009l iMeneghetti et al.||201Ua| |201ip. A handful 



of the clusters in the CLASH A-ray selected subset have 
some evidence for departures from symmetric X-ray sur- 
face brightne ss d istributions. These systems are briefly 
discussed in §3.3[ 

These 20 clusters comprise the primary sample for our 
DM distribution studies. Existing X-ray data will be 
used to measure the gas density profile fo r subtraction 
from the total mass profile (as applied by Lemze et al. 
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(2008) for A1689), which is necessary when making de- 
t ailed comparison with the predictions of DM-only sim- 
ulations. Although these clusters are X-ray selected, one 
or more lensed arcs are indeed visible in all of the clusters 
for which sufficiently high-quality imaging was available. 
This indicates the X-ray selected clusters in our sample 
have Einstein radii in the range 15" to 35". X-ray im- 
ages from the ACCEPTj^ database for these 20 CLASH 
clusters are shown in Figure [6] 

In addition, CLASH includes five clusters selected 
based solely on their exceptional strength as gravita- 
tional lenses (large Einstein radii, Oemi > 35" for = 2). 
Analysis of these clusters will allow us to further quan- 
tify the lensing selection bias toward high concentrations. 
The excellent strong lensing data yielded by these clus- 
ters will enable derivation of some of the highest re solu- 
tion d ark matter maps which may be obtained (as in |Coe| 
et al.| 2010). The primary motivation for selecting these 
five "nigh-magnification" clusters, however, was to signif- 
icantly increase the likelihood of discovering very highly 
magnified high redshift galaxies (i ]2.3[ ). While most of 
the CLASH clusters lie behind regions of low Galactic 
extinction (median E(B-V) = 0.026), there is one cluster, 
MACS2129. 4-0741, with significant amounts of Galactic 
cirrus located northeast of the cluster core. Because this 
cluster has one of the larger Einstein radii known, its in- 
clusion in the "high-magnification" sample was deemed 
worthwhile in spite of the presence of this cirrus. 

Figure [7] shows, respectively, the distribution of the 
projected separation between the BCG and the peak of 
the X-ray surface brightness and the distribution of the 
ellipticity of the X-ray emission f or the CLASH sample . 
EUipticity measurements are from Maughan et al. ( 2008 ). 
The X-ray selected cluster sub-sample exhibits, on aver- 
age, a very small off'set between the BCG and peak X-ray 
flux. The X-ray selected cluster sub-sample also exhibits 
a lower average ellipticity than our high-magnification 
(Einstein-radii selected) sample but is consistent with 
the mean ellipticity of the larger Maughan et al. ( 2008 1 
sample. 

The 25 clusters in the CLASH program are presented 
in Table [3] and some of their key X-ray properties are 
given in Table |4j including the bolometric luminosity 
(defined for convenience to be between 0.1 and 100 keV), 
temperature, and cluster-to-solar [F e/H] ratio, where the 
solar abundance reference is from lAnders fc Grevessel 
( [1989) . Table |4[ also lists the source of the X-ray selec- 
tion m column D. The CLASH sample is dr awn heavily 



from the AbeU and MACS cluster catalogs ( Abefl 1958 
Abefl et al.|[T989l [Ebehng et al~1[200T| [20071 [MU| ^TEe 



CLASH clusters span almost an order ot magnitude in 
mass (- 5 to - 30 X lO^^^ Mq). These clusters were also 
selected to cover a wide redshift range (0.18 < z < 0.90 
with a median z^ed — 0.40) allowing us to probe the full 
c(M, z) relations expected from simulations. 

The X-ray parameters in Table [4| are derived by us 
using CIAO v4.3 and CALDB v4.4.3. We filtered the 
Chandra data for fiares and reprojected the deep Chan- 
dra background fields to match the cluster observations. 
The background data were also filtered for "status^ 0" 

25 ACCEPT = Archive of Chandra Cluster Entropy Profile Ta- 
bles. See http:/ /www. pa. msu.edu/astro/MC2/accept/ for more in- 
formation 



events to be suitable for use with the VFAINT mode data. 
To obtain X-ray luminosities and temperatures uniformly 
across the sample, we extracted spectra from apertures 
with radii of 714 kpc (500 H^qq kpc) and excluded the 
central 71.4 kpc (50 /i|j~oo kpc). These "core-excised" 
spectra provide X-ray temperature estimates that are rel- 
atively u naffected by the pr esence or absence of a cool 



core (e.g., Markevitch ( 1998[ )). Figure [8| shows the X-ray 
temperature, l'^;, as a function of the scaled bolomet- 



ric luminosity, Lj,/i?(z), for CLASH clusters using the 
data in Table ^ The X-ray luminosities are scaled as- 
suming self-similar evolution, E{z) — ^^^^(l + z)^ + A, 
to allow direct comparison to th e low-redshift luminosity 
temperature relat ionships from [Markevitch ( [l998 ) and 
Pratt et al l ( |2009[ ) (La (measured within 0.15-XDD R500) 
vs. 'i'3 (0.15-0.75 R500), to be specific.) We verified that 
the deep background particle event rates matched that 
seen in the cluster spectra between 9—12 keV. The X- 
ray counts were binned to a minimum of 20 counts per 
energy bin. Since the extracted spectra were dominated 
by source counts, the results were not very sensitive to 
the background scaling. To estimate the cluster temper- 
ature, we used XSPEC 12.6.0q to fit the X-ray spectra 
between 0.7 — 7.0 keV. We assumed a single tempera- 
ture plasma model (the XSPEC apec model) absorbed by 
a Galactic hydrogen c olumn fixed to the value ob tained 
from the Bell survey ( [Dickey fc Lockmanj [l990| ) . The 
metallcity, temperature, and normahzation were allowed 
to be free. The best fit temperatures and la error bars, 
and the bolometric X-ray luminosities are reported in 
Table H 

3.2. Cluster Sample Size Requirements 

The required size of our "relaxed" cluster sample is de- 
rived from the goal to measure "average" cluster concen- 
trations to ^ 10% (after accounting for variations in mass 
and redshift) and to detect a ^ 15% deviation from the 
concentrations of simulated clusters at 99% confidence. 

Statistically, if we assume measurement errors are in- 
dependent and normally distributed then the number of 
clusters, N^usj needed to measure the average concen- 
tration to a fractional accuracy of / is just 



(1) 



where the total scatter in an individual measurement of 
the concentration, atot, is 



2 

^tot 



2 

"int 



■ (2) 

neas \ J 

and includes contributions from intrinsic variation (in- 
cluding the effects of orientation averaging of triaxial 
halo shapes), variations due to intervening large scale 
structure (LSS), and measurement uncertainties. Con- 
centrations of relaxed simu lated cluster halos have intrin- 
sic scatters of ^ 30% (e. g., Neto et al.|2007 Maccio et al. 
2008 Duffy et al. 2008). LSS can scatter conc entration 
measurements by 13% ( Hoekstra et al.|2002 ). And we 
conservatively estimate measurement uncertainties to be 
~ 37% as follows based on empirical data (though analy- 
sis of simulated data sug gests ^11% accuracy is possible 
Meneghetti et al.||2010a[ ). 

Analysis ot CLij024-|-17 yielded a concentration mea - 
surement uncertainty of ~ 22% (Umetsu et al. 2010), 
which we will use here as a baseline. We assume here 
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Figure 6. Cutouts of Chandra X-ray images centered on the 20 CLASH clusters in the X-ray selected subsample. These images are 
taken from the Archive of Chandra Cluster Entropy Profile Tables (ACCEPT). Each cutout subtends ~3.45 arcminutes, which is nearly 
the same as the ACS/WFC field of view on HST. 

that our precisions will be dominated by the quality 
of our SL data and that our measurement uncertainty 
will scale roughly as cr cx 1/V Narc c< 1/Re, where for 
CL0024-I-17 with i?B - 35", Narc = 33 arcs were de- 
tected. Note we have also assumed Narc oc R%. Based on 
archival WFPC2 and ACS images, we measured a mean 
Re ^ 20" for our sample of 20 relaxed clusters, which 
should conservatively yield Narc ^ 12 arcs per cluster. 
This scales the ~ 22% uncertainty for CL0024-I-17 with 
33 arcs up to a conservative ^ 37% for an average cluster 
with just 12 arcs. 

Adopting (Tint = 0.30, aLSS = 0.13 and amcas = 0.37, 
we find atot = 0.49, and hence, if / = 0.10 then 



Nr., 



24. This empirical sample size estimate is con- 



sistent with one derived from numerical simulations of 
strong lensing, as shown in Figure |9] The approximately 
log-normal distributi on of DM halo concentra tions seen 
in such simulations ( Meneghetti et ah] 2010b) indicates 
that ~20 massive clusters (M^,^ > 5 x lO^'' M©) are re- 
quired to detect 15% deviations from expectations in the 
standard cosmological model at 99% confidence. 

Note that we may expect deviations to be greater than 
15% if results from lensing-selected clusters are any indi- 
cation. Simulated IO^^Mq/i"^ clusters at z = typically 



have Cyi r 5 (e.g., lDuffy et al. 20081. For this fiducial 
cluster, |Ogur i et al.^ (2009) instead find Cyir ^ 12 for 
observed clusters (based on a fit to values observed for 
10 lensing clusters accounting for the various masses and 
redshifts). Even after including a 50% lensing bias, these 
observed concentrations are still ^ 70% greater than ex- 
pectations. Being more conser vative and assuming a fac - 
tor of 2 (100%) lensing bias fMeneghett i et al.| |2010a), 
the observed concentrations would still be ^ 2 0% g reater 
than expectations. Although as disc ussed in |2.1[ more 
recent simulations (Prada et al. 2011 ) may alleviate these 
discrepancies. 

3.3. Notes on Clusters with Possible Substructure 

While our X-ray selection criteria favors the inclusion 
of highly relaxed clusters in the CLASH sample, for 8 of 
our clusters the dynamical state is somewhat ambiguous. 
Some researchers have reported evidence for substructure 
in the X-ray surface brightness profiles of these 8 clusters. 
The presence of substructure may suggest that a cluster 
is not fully dynamically relaxed. We summarize this evi- 
dence below. We note that the presence of a minor degree 
of substructure does not inhibit our ability to determine 
a cluster's mass profile characteristics. In simulations. 
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Figure 7. Left: The distribution of the projected separation between the Brightest Cluster Galaxy (BCG) and the peak of the X-ray 
surface brightness for the X-ray selected (dark grey) and Einstein-radii selected (red) cluster subsamples. The average projected BCG - 
X-ray peak separation for the X-ray selected sample is 10.4 kpc. Right: The distribution of the ellipticity of the X-ray emitting intracluster 
gas for the X-ray sele cted (dar k grey) and Einstein-radii selected (red) cluster subsamples. The distribution of the ellipticities of the 96 
clusters in the Maugh an et al.| ||2008 ) sample that are not in common with CLASH is shown in the light grey histogram. The mean X-ray 
emission ellipticities tor tue A-ray and Einstein-radii selected cluster samples are 0.19 and 0.28, respectively. 
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Abell 209 


01:31:52.57 


-13:36:38.8 


0.206 


19 


20 


12451 


(3) 


Abell 383 


02:48:03.36 


-03:31:44.7 


0.187 


18 


20 


12065 


(3) 


MACS0329.7-0211 


03:29:41.68 


-02:11:47.7 


0.450 


19 


20 


12452 


(2.5) 


MACS0429.6-0253 


04:29:36.10 


-02:53:08.0 


0.399 


20 


20 


12788 


(0.5) 


MACS0744.9-f3927 


07:44:52.80 


+39:27:24.4 


0.686 


18 


17 


12067 


6 


Abell 611 


08:00:56.83 


+36:03:24.1 


0.288 


19 


18 


12460 


2 


MACS1115.9-f0129 


11:15:52.05 


+01:29:56.6 


0.352 


19 


20 


12453 


0.5 


Abell 1423 


11:57:17.26 


+33:36:37.4 


0.213 


20 


20 


12787 


(0.5) 


MACS1206.2-0847 


12:06:12.28 


-08:48:02.4 


0.440 


18 


20 


12069 


0.5 + (0.5) 


CLJ1226.9-I-3332 


12:26:58.37 


+33:32:47.4 


0.890 


20 


18 


12791 


16 


MACS1311.0-0310 


13:11:01.67 


-03:10:39.5 


0.494 


20 


20 


12789 





RXJ1347.5-1145 


13:47:30.59 


-11:45:10.1 


0.451 


18 


15 


12104 


6 + (0.5) 


MACS1423.8+2404 


14:23:47.76 


+24:04:40.5 


0.545 


20 


17 


12790 


5 


RXJ1532.9-f3021 


15:32:53.78 


+30:20:58.7 


0.345 


19 


20 


12454 


(0.5) 


MACS1720.3-f3536 


17:20:16.95 


+35:36:23.6 


0.391 


19 


20 


12455 


(0.5) 


Abell 2261 


17:22:27.25 


+32:07:58.6 


0.224 


18 


20 


12066 


(0.5) 


MACS1931.8-2635 


19:31:49.66 


-26:34:34.0 


0.352 


19 


20 


12456 





RXJ2129. 7+0005 


21:29:39.94 


+00:05:18.8 


0.234 


19 


20 


12457 


(1) 


MS2137-2353 


21:40:15.18 


-23:39:40.7 


0.313 


18 


18 


12102 


5 + (7) 


RXJ2248. 7-4431 (Abell 1063S) 


22:48:44.29 


-44:31:48.4 


0.348 


19 


20 


12458 


(0.5) 


High Magnification Clusters; 
















MACS0416. 1-2403 


04:16:09.39 


-24:04:03.9 


0.42'' 


19 


20 


12459 


(1) 


MACS0647.8+7015 


06:47:50.03 


+70:14:49.7 


0.584 


18 


18 


12101 


9 


MACS0717.5+3745 


07:17:31.65 


+37:45:18.5 


0.548 


18 


17 


12103 


7 


MACS1149.6+2223 


11:49:35.86 


+22:23:55.0 


0.544 


18 


18 


12068 


5 


MACS2129.4-0741 


21:29:26.06" 


-07:41:28.8" 


0.570 


18 


18 


12100 


5 



^ Central cluster coordinates derived from optical image instead of X-ray image 

^ Cluster redshift for MACS0416. 1 — 2403 is based on Chandra X-ray spectrum (this work). Uncertainty on this value is 
±0.02. For this fit, the core was not excluded. To maximize the counts, an aperture of 714 kpc was used, binned to achieve 
a minimum of 20 counts per energy bin. The background in the 0.7-7.0 keV range was fit from the Chandra deep fields. 

Archival ACS or WFPC2 imaging data only; WFPC2 orbits shown in parentheses. Archival ACS images, when available, 
are used in conjunction with new CLASH data to achieve the desired depths in all filters. 
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Table 4 

X-ray Properties of the CLASH Cluster Sample 





kT 




[Fe/H] Ratio 


Ellipticity'' 


Centroid Shift'' 


X-ray Morphology 


Cluster 


(keV) 


(10*-* erg s-i) 


(Cluster/Solar)<= 


e 




Reference 


X-ray Selected Clusters: 














Abell 209 


7.3 ± 0.54 


12.7 ±0.3 


0.18 ±0.09 


0.21 ± 0.01 


0.55 ±0.05 


1 


Abell 383 


6.5 ± 0.24 


6.7 ±0.2 


0.59 ±0.10 


0.04 ± 0.01 


0.18 ±0.02 


2 


MACS0329. 7-0211 


8.0 ± 0.50 


17.0 ±0.6 


0.48 ±0.10 


0.15 ±0.03 


1.38 ±0.13 


2,3 


MACS0429.6-0253 


6.0 ± 0.44 


11.2 ±0.5 


0.34 ±0.11 


0.21 ± 0.02 


0.38 ±0.03 


2,3 


MACS0744.9-I-3927 


8.9 ± 0.80 


29.1 ± 1.2 


0.37 ±0.10 


0.11 ± 0.02 


1.41 ± 0.13 


2,3 


Abell 611 


7.9 ± 0.35 


11.7 ±0.2 


0.30 ± 0.07 






2 


MACS1115.9-I-0129 


8.0 ± 0.40 


21.1 ± 0.4 


0.32 ± 0.06 






2,3 


Abell 1423 


7.1 ± 0.65 


7.8 ±0.2 


0.23 ±0.13 






4 


MACS1206. 2-0847 


10.8 ± 0.60 


43.0 ±1.0 


0.25 ± 0.09 






3,4 


CL J 1226. 9-1-3332 


13.8 ± 2.80 


34.4 ± 3.0 


0.36 ±0.25 


0.10 ± 0.03 


1.30 ±0.11 


2 


MACS1311.0-0310 


5.9 ± 0.40 


9.4 ±0.4 


0.42 ± 0.10 


0.08 ±0.02 


0.39 ±0.03 


2,3 


RXJ1347.5-1145 


15.5 ± 0.60 


90.8 ±1.0 


0.20 ±0.06 


0.26 ±0.01 


0.63 ±0.05 


2 


MACS1423.8-f2404 


6.5 ± 0.24 


14.5 ±0.4 


0.35 ± 0.06 


0.17 ± 0.03 


0.25 ±0.02 


2,3 


RXJ1532. 9-1- 3021 


5.5 ± 0.40 


20.5 ± 0.9 


0.52 ± 0.12 


0.18 ±0.02 


0.07 ± 0.01 


2 


MACS1720.3-I-3536 


6.6 ± 0.40 


13.3 ±0.5 


0.29 ±0.09 


0.17 ± 0.02 


0.57 ± 0.05 


2,3 


Abell 2261 


7.6 ± 0.30 


18.0 ±0.2 


0.31 ± 0.06 


0.10 ± 0.01 


0.71 ± 0.06 


5 


MACS1931.8-2635 


6.7 ± 0.40 


20.9 ± 0.6 


0.26 ±0.08 


0.30 ± 0.01 


0.28 ±0.02 


2,3 


RXJ2129. 7-1-0005 


5.8 ± 0.40 


11.4 ±2.0 


0.33 ±0.10 


0.26 ±0.02 


0.55 ±0.05 


2 


MS2137-2353 


5.9 ± 0.30 


9.9 ±0.3 


0.41 ± 0.08 






2 


RXJ2248.7-4431 (Abell 1063S) 


12.4 ± 0.60 


69.5 ± 0.1 


0.39 ± 0.06 


0.20 ± 0.01 


0.74 ± 0.06 


6 


High Magnification Clusters: 














MACS0416. 1-2403 


7.5 ± 0.80 


16.0 ±0.9 


0.40 ±0.14 






3 


MACS0647.8-I-7015 


13.3 ± 1.80 


32.5 ± 2.1 


0.30 ±0.19 


0.36 ±0.02 


0.62 ±0.06 


3 


MACS0717.5-I-3745 


12.5 ± 0.70 


55.8 ± 1.1 


0.22 ± 0.07 


0.30 ± 0.01 


2.11 ± 0.18 


3 


MACS1149.6-I-2223 


8.7± 0.90 


30.2 ± 1.2 


0.24 ±0.11 


0.25 ±0.02 


1.21 ± 0.10 


3 


MACS2129.4-0741 


9.0 ± 1.20 


22.6 ± 1.5 


0.40 ±0.17 


0.19 ±0.03 


1.56 ±0.14 


3 



Ref erences. — (1) |Maughan et al.| | |2008t ; (2) [Alien et al.| 1 (2008} ; (3) [Ebeling et al.||2007^ ; (4) [Cavagnolo et al.] l |2008| l; (5) |Mantz et al.| 
< |2010( ; (6) CLASH team selected cluster. 

^ The X-ray bolomctric luminosity covers the energy range 0.1 — 100 keV. 

^ The X-ray cllipticity, e, and centroid shift, (w) , values a re from [MaAJ^h^Ji^^taL ||2008| . 
^ The solar abundance reference values used here are from [Anders iS^Grevesse |l |1989[ l. 
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Figure 8. The X-ray temperature, Tx , as a function of the scaled 
bolometric luminosity, /E(z), for the CLASH cluster sample. 
Data points are from Table |4] The X-ray luminosities are scaled 
assuming self-similar evolution, to allow direct comparison to the 
low-rcdshif t luminosity tempe rature relationships from Markevi tchI 
JTqQSJ andjPratt et aL]l|2009|). Th e - relations from Marg] 
|viCT[l|1998p and lPratt erar[ T |2009[ l (L2 (measured within 0.15-l.UU 
tl 50oj VS. T3 (U.i5-U.Y5 Ksoojj are shown. CLASH clusters are 
representative of these larger X-ray cluster samples. 




number of clusters 

Figure 9. Mean "over-concentration" ratio necessary to reject 
with 99% confidence the hypothesis that measured concentrations 
are drawn from a sample with the expected values. Average de- 
viations of ^ 15% from expectations will be detected with 99% 
confidence based on our relaxed sample of 20 clusters. We have 
adopted the finding from ACDM simulations that DM halo con- 
centrations are log-normal ly distributed with a standard deviation 
of (T Si 0.25 ± 0.03 (e.g., [Duffy et ari|2008| . The shaded band 
includes this uncertainty. 

relaxed clusters have DM concentrations that are, on av- 
erage, ^ 20% higher th an the g eneral cluster population 



(Duffy et al.[2008.; Prada et al. 



2011). 



Abell 209: This system is ma rginally unrelaxed ac 



cording to Smith et al. ( 2005 ) . They based their re- 
laxed/ unr5a5ce3"cEajacterIzat ion mainly on a mass 
ratio between the cluster main component and the 
total mass (both at r < 500 kpc). Clusters with 
mass ratio values below 0.95 are considered unre- 



Their value for Abell 209 is 0.87 ± 0.06 
hereafter 



Maughan et al 
an elhpticity o f e 



(|2008|, hereafter M08, derived 
^■1^21 ± 0.01 for this cluster. 



Gilmour et al.| ( |2009 ), hereafter G09, find it to be 
relaxed based on a visual examination of its X-ray 
morphology. 



• MACSJ0329.7-0211 and CLJ1226.9-H3332: these 
two clust ers are included in t he A0 8 compilation (as 
well as in [Schmidt fc Allen| ( |2007l ) , hereafter SA07, 
and were classihed in these works as dynamically 
relaxed. However, M08 report that both exhibit 
evidence for substructure. 



• MACSJ0744.9-I-3927: This cluster shows evidence 
for substructure (SA07; M08). On the other hand 
G09 find it to be relaxed based on a visual exami- 
nation of its X-ray morphology. 

• MACSJ1206.2-0848: In both optical and X-ray im- 
ages the cluster appears close to relaxed in projec- 
tion, with a pronounced X-ray peak at the location 
of the BCG. However, some evidence of merger ac- 
tivity along the line of sight may be suggested by 
the very high velocity dispersion of 1580 km s~^. 
G09 visually classify this cluster as relaxed. 

• RXJ1347. 5-1145: Significant substructure is ob- 
served in this system (SA07; M08). M08 mea- 
sure an average elhpticity of the X-ray isophotes 
of e = 0.26 ±0.01. 

• Abell 2261: M08 found a small level of substructure 
and G09 classified it as disturbed. 

• RXJ2248. 7-4431: M08 found this cluster to be 
slightly elliptical e = 0.2 ±0.01 and with some level 
of substructure, while G09 found it to be relaxed. 

4. SURVEY DESIGN AND IMPLEMENTATION 

The CLASH program consists of 524 HST orbits, in- 
cluding 50 for supernova follow-up. The bulk of the 
program (474 orbits) will be used to image 25 galaxy 
clusters each to a depth of 20 orbits divided among 16 
HST /ACS and WFC3 filters as shown in Table [5| and 
Figure [TOj The number of required cluster orbits is re- 
duced from 500 to 474, as some of these filter depths 
have already been achieved in existing data. Columns 6, 
7, and 8 in Table [3] provide the number of new orbits al- 
located, the 77S'T Program ID (where available), and the 
number of existing archival orbits available, respectively, 
for each cluster. The multi-band observations span the 
near-ultraviolet to near- infrared (2,000 - 17,000A). In 
the wavelength range ^4, 000-9,000A, we use A CS for its 
greater throughput efficiency (especially toward the red 
end) and larger observing area. 

A typical CLASH observing sequence is presented in 
Tabled Each cluster is observed at two orientations to 
increase the area covered in the parallel field SN search. 
The two orientation angles are typically ^ 30° apart to 
minimize the overlap between the ACS parallel point- 
ings 



(see 



4.3 for more details). Although we label the 



visits as "A" and "B" in Table [6j either orientation may 
be executed first. The choice depends solely on schedul- 
ing constraints. Indeed, there are often overlaps in time 
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Figure 10. Each CLASH cluster is observed in 16 HSr filters spanning ~2,000-17,000A with WFC3/UVIS in the near-ultraviolet, ACS 
in the optical (extending into the near-infrared), and WFC3/IR in the near-infrared (see Table[5]l. Total throughput curves are plotted for 
each filter. For clarity, some curves are offset vertically by 0.2 (dashed line). 



Table 5 

CLASH Exposure Times, Limiting Magnitudes, and Extinction Coefficients 



Camera / 
Channel 


Filter 
Element 


Orbits 


Average Exposure 
Time (scc)° 


lOcr Limit 

(AB mag)*" 


5o" Limit 
(AB mag)*" 


Galactic Ext. 
(AB mag/E(B-V)) 


WFC3/UVIS 


F225W 


1.5 


3558 


25.7 


26.4 


7.474 


WFC3/UVIS 


F275W 


1.5 


3653 


25.7 


26.5 


6.140 


WFC3/UVIS 


F336W 


1.0 


2348 


25.9 


26.6 


5.090 


WFC3/UVIS 


F390W 


1.0 


2350 


26.5 


27.2 


4.514 


ACS/WFC 


F435W 


1.0 


1984 


26.4 


27.2 


4.117 


ACS/WFC 


F475W 


1.0 


1994 


26.8 


27.6 


3.724 


ACS/WFC 


F606W 


1.0 


1975 


26.9 


27.6 


2.929 


ACS/WFC 


F625W 


1.0 


2008 


26.4 


27.2 


2.671 


ACS/WFC 


F775W 


1.0 


2022 


26.2 


27.0 


2.018 


ACS/WFC 


F814W 


2.0 


4103 


27.0 


27.7 


1.822 


ACS/WFC 


F850LP 


2.0 


4045 


26.0 


26.7 


1.473 


WFC3/IR 


F105W 


1.0 


2645 


26.6 


27.3 


1.015 


WFC3/IR 


FllOW 


1.0 


2415 


27.0 


27.8 


0.876 


WFC3/IR 


F125W 


1.0 


2425 


26.5 


27.2 


0.757 


WFC3/IR 


F140W 


1.0 


2342 


26.7 


27.4 


0.609 


WFC3/IR 


F160W 


2.0 


4920 


26.7 


27.5 


0.470 




WFC3: UVIS IR 



^ Exposure times are the average time for each filter for all CLASH cycle 18 observations. Exposure times for 

cycles 19 and 20 may differ slightly due to scheduling considerations. 

^ Limiting magnitudes are for a circular aperture that is 0.4 arcsec in diameter. 

Th e su rvey footprint and exposure map arc shown in Fig- 
ure [TT] Color images of the co-added H ST d&ta. from the 
first 4 CLASH cluster are shown in Figure [T2j 

CLASH clusters will be observed over tHe course of 
three annual HST observing cycles, with 10, 10, and 5 
clusters to be observed in cycles 18, 19, and 20, respec- 
tively. The current assignments of the cluster observa- 
tions to the HST cycle number are shown in column 5 of 
Table [3j The duration of this MCT program is driven, in 
part , by the cadence required by supernova search (see 
S}43] below). 

4.1. 




ACS 



Figure 11. Each cluster is observed at two orientations approx- 
imately 30 degrees apart to increase the supernova search area in 
the parallel fields. WFC3 parallel observations arc obtained while 
the cluster core is being observed with ACS and vice versa. At 
the median cluster redshift of Zyj^^^ = 0.4, the 6 arcminute separa- 
tion between the center of the parallel field and the cluster center 
corresponds to a projected distance of 1.9 Mpc. 

when the "A" and "B" orientations are both being ex- 
ecuted. When the entire sequence of exposures for a 
cluster is completed, the region covered by all 16 filters 
subtends an area of 4.08 square arcminutes. Larger areas 
about the cluster center are covered by ACS in 7 filters. 



Filter Selection and Exposure Times 

Redshift estimates for multiply-lensed images are cru 
cial for breaking lensing degener acies and tightening con 
straints on mass profiles (e.g ., iBroadhurst et al.||2005 



Zitrin et al 
ot the usetu 



2009 

lensed images are much' 



Saha & Read 20091). However, most 
h too faint for spec- 



troscopy. Typical lensed source magnitudes are 23 < / < 
28 (Figure 131, so that only the brightest arcs yield spec- 
troscopic redshifts even when observed with the largest 
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Table 6 

Typical CLASH Observing Sequence 







Primary 


Primary Cluster 


Parallel 


Parallel Field 


Visit ID 


Epoch'' 


Camera 


Target Filters 


Camera 


Filters 


Orientation A: 












AO, Al 


lA 


ACS 


F625W, F850LP 


WFC3 


F350LP, F125W, F160W 






F475W, F775W 






A2 


lA 


WFC3 


FllOW, F160W 


ACS 


F775W, F850LP 


A3 


2A 


ACS 


F606W, F814W 


WFC3 


F350LP, F125W. F160W 


A4 


2A 


WFC3 


FllOW, F160W 


ACS 


F775W, F850LP 


A5 


3A 


ACS 


F435W, F814W 


WFC3 


F350LP, F125W, F160W 


A6" 


3A 


WFC3 


F225W, F390W 


ACS 


F775W, F850LP 


AO 




A 






r oOUij-r, r IZO VV . -T IDU VV 


A9 


4A 


WFC3 


FllOW, F160W 


ACS 


F775W, F850LP 


Orientation B: 












BO 


IB 


ACS 


F625W, F850LP 


WFC3 


F350LP, F125W, F160W 


HI 


IB 


WFC3 


F125W, F160W 


ACS 


F775W, F850LP 


B2 


2B 


ACS 


F606W, F775W 


WFC3 


F350LP, F125W, F160W 


B3" 


2B 


WFC3 


F275W, F390W 


ACS 


F775W, F850LP 


B5 


3B 


ACS 


F435W, F814W 


WFC3 


F350LP, F125W. F160W 


B6 


3B 


WFC3 


F336W, F105W, F140W 


ACS 


F775W, F850LP 


B7 


4B 


ACS 


F475W, F850LP 


WFC3 


F350LP, F125W, F160W 


B8 


4B 


WFC3 


F125W, F160W 


ACS 


F775W, F850LP 


Visits A6 and B3 


arc each executed as 


a single 2-orbit visit using 


the same g 


uidc star acquisition. This is 



essential to facilitate alignment of the NUV exposures. 

^ Epochs lA — 4A are sequential in time, as are epochs IB — 4B. However, depending on scheduling constraints, 
exposures for orientation B may precede or be interlaced with those for orientation A. 



ground-based facilities. With continuous sampling of the 
broad wavelength range from the NUV to NIR (~ 2,000- 
17,000 A) that is enabled with WFC3 and ACS we can 
now obtain very accurate photometric redshifts (photo- 
z's) for most of the lensed objects down to an apparent 
F775W AB magnitude limit of 26. 

For a fixed total observing time, splitting observations 
into multiple (ideally overlappin g) filters significantly im- 
proves photo-z precision ( ,Benftez et al. 2009) . We per- 
formed simulations to inform our hlter selection and 
estimate our eventual photo-z precision. Galaxy mag- 
nitudes, redshifts, and SEDs (spe ctral energy dist ribu- 
tions) were drawn from the UDF (Coe et al. 2006) and 

adding noise as appro- 
Photo- 



"re-observed" with our filter set 
priate given our proposed depths in each filter 
z's were then re-estimated using BP Z (Bem'tez 



Benftez et al 
tion, we hnd 



[20041 |Coe et al.||2006| . 

hat for 16 hlters, we obt ai 



2000 
In this Simula- 



we need to achieve a lOcr limiting AB magnitude of 26 
in each filter. We augment this by two orbits per clus- 
ter to extend the NIR depth to a lOcr limiting AB mag 
of F160W = 26.7 and F814W = 27.0 as needed for our 
lensed high-z galaxy search (j |2^3| and see Tabled. The 
limiting magnitudes in Table are for a 0.4" diameter 
circular aperture and a point source with a flat J^^ spec- 
trum. The minimum exposure in any given filter is one 
orbit and observations in each filter are distributed over 
at least four exposures to ensure robust cosmic ray rejec- 
tion. 

The five NIR filters provide the ability to identify 
z > 7 galaxies with high confidence. Measured colors in 
these filters are essential for discriminating between such 
high-redshift galaxies and lower-redshift highly-reddened 
galaxies (e.g., see Figure 14). The former are generally 



am very accurate 
{Az ~ 0.02(1 + z)) photo-z's for 80% of objects with 
F775W mag < 26. Most importantly, we find that we 
will be able to acquire ~ 6 times as many reliable photo- 
metric redshifts than spectroscopic redshifts for objects 
at z > 1, enabling a very substantial improvement in the 
number of uniqu e co nstraints on the DM mass distribu- 
tions (see Figure 13). 

The coverage provided by the 16 WFC3 and ^C5fihers 
allows the Lyman-limit feature (rest frame 912A) to be 
photometrically traced to redshifts as low as z ~ 1.5 and 
Lya to be detected out to z ^ 10. The inclusion of 
NUV photometry, for example, resolves one of the most 
common photo-z degeneracies between the Balmer break 
in z ^ 0.2 galaxies an d the Lyman break in z ~ 3 galaxies 
(IRafelski et al. 112009 1. 



relatively flat in (or AB magnitudes) in the NIR while 
the latter generally increase in brightness as a function of 
wavelength. Quantitatively, the following selection crite- 
ria have proven effective: z — J > 0.8 (a flux decrement 
factor of 2.1) 



et al. 2009 



and J-H < 0.5 (Bradley et al.|2008,; Zheiigl 


Bradley et al. 


2011 


). lb reliably measure 



'I'he exposure times for the primary camera (cluster 
center position) are set primarily by the need to achieve 
80% photometric redshift completeness down to F775W 
= 26 AB mag. This requires ^18 orbits per cluster as 



the above colors for high-z lensed galaxy detection, our 
exposures are set to reach the following AB magnitude 
limits: F850LP: 26.7 (5ct), FllOW: 27.8 (5cr), F125W: 
26.5 (lOcr), and F160W 26.7 (lOcr). 

4.2. Dither Pattern 

In each orbit we use a compact 4— point dither pat- 
tern that provides half-pixel sampling along both detec- 
tor axes. The dither pattern serves to both improve the 
spatial sampling of the point spread function (PSF), es- 
pecially for the WFC3/IR detector with its large pixel 
scale of 0.128 arcsec/pixel, and to help remove hot pix- 
els and other detector imperfections that may be un- 
accounted for in the calibration reference files. In the 
first epoch of each roll angle, we use a small-scale dither 
pattern to preserve, as best as possible in light of signifi- 
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Figure 12. Cutouts of the CLASH images of tlie central regions of (A) Abell 383, (B) MACS1149.6+2223, (C) Abell 2261, and (D) 
MACS1206.2— 0847. The field of view shown in each image is 65 X 65 arcseconds and is centered on the brightest cluster galaxy. These 
color composites are made from a combination of 12 of the 16 filters available for each cluster. The blue channel is a sum of the ACS F435W 
and F475W filters, the green channel is a sum of the ACS F606W, F625W, F775W, F814W, and F850LP filters, and the red channel is a 
sum of all five of the WFC3/IR broad band filters. Each cutout shown here represents ^ 28% of the area covered by all 16 filters and just 
~ 10% of the total area imaged in the cluster core. 

Parallel observations are being obtained with a pri- 
mary science goal of detecting Type la SNe. While the 
cluster core is being imaged with ACS, we obtain par- 
allels with WFC3, and vice versa. For each of the two 
roll angles, observations are obtained over the course of 
four epochs spread out over ^ 30 to 45 days 10 to 
14 days between observations). Image differencing anal- 
ysis enables us to detect supernovae which have gone off 
between observations. As noted above, the exposure se- 
quences are executed either sequentially or interleaved, 
depending on scheduling availability. 

Upon detection of a promising high-z type la SN in 
the parallel field around a cluster, we can reprogram the 
remaining parallels to provide the follow-up (light curve 



cant geometric distortions, half-pixel sampling across the 
detector. In subsequent epochs involving WFC3/IR ob- 
servations, either in prime or parallel, we use a slightly 
larger dither pattern to help identify and remove per- 
sistence artifacts from compact sources, which, if uncor- 
rected, could possibly be misidentificd as SNe candidates. 
While our small-scale dither patterns are much smaller 
than is needed to cover the WFC3/UVISimA ACS/WFC 
detector gaps, our cluster observations are obtained at 
two orientations, leaving only two small diamond-shaped 
regions (~ 4.4 square arcseconds each) in the central clus- 
ter area without data in all 16-filters. 

4.3. Observation Cadence and Supernova Follow-Up 
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Figure 13. Top: Magnitude distribution of 132 multiply lensed 
images detected in A1689 and CL0024+17 fBroadhurst et al.pOOS] 
[Zitrin et al. 2009). Most are too faint for spectroscopic lollow-up! 
ubttom: tiy increasing the number of filters for a fixed observing 
time, we show how photometric rcdshifts are improved using two 
metrics. The horizontal axis shows, for example, that 16 filters 
yield very accurate ( ~ 0.02(1 + z)) photo-z's for 80% of F775W 
mag < 26 objects. This allows us to obtain precise redshifts for ~ 6 
times more objects than attainable using spectroscopy, as shown 
on the vertical axis. These estimates are based on simulated pho- 
tometric and photo-z cata logs derived from galaxies observed in 
the UDF | |Coe et al.|2006t . 



imaging and grism spectroscopy) necessary to measure 
its distance. The ability to reprogram the later of the 
two orientations to follow-up a SN detected early in a 
cluster observing sequence sets the upper limit on the 
angular offset between the two orientations - both ori- 
entations must be accessible during the entire cluster se- 
quence. This constraint means that the two orientations 
cannot be more than ~ 30° apart. A lower limit on the 
angular separation between the two orientations is set 
to ^ 20° to ensure there is not excessive overlap in the 
ACS parallel pointings (with the exception of RXJ1347 
where an orientation shift of 15° was needed to avoid 
placing a bright star in the parallel fields). The abil- 
ity to reprogram a downstream orientation shift means 
that a new target of opportunity (ToO) observation will 
not be required if the SN is discovered early in the clus- 
ter observing sequence. In this regard, the design of the 
CLASH SN program has part of the follow-up built into 
its implementation. For the flexibility required to fol- 
low SNe found at the end of their corresponding cluster 
observing sequence, a small number of ToOs (four) and 
follow-up reserve orbits are included in the CLASH or- 
bit allocation. CLASH parallel fields are not observed in 
16 filters. The ACS parallels are taken in F775W and 
F850LP and the WFC3 parallels are taken in F350LP 

(UVIS), F125W (IR) and F160W (IR) (see Table §). 

The CLASH and C ANDELS (Grogin et aTj^^oTT 



200 orbits: 50 from the CLASH allocation and 150 from 
CANDELS. The coordinated program is led by A. Riess 
who is a co-I in both programs. 

4.4. ACS Failure Options 

While ACS functionality was restored in SM4, it is 
now only a "single-string" instrument, meaning there is 
no redundant path if the main CCD electronics box ex- 
periences a failure. This was a constraint imposed by 
the nature of the ACS repair. The primary impact of 
such a failure would be the loss of our parallel obser- 
vations for supernova searches. Supernovae could still 
be detected in our primary cluster core observations as- 
suming we were to continue to distribute the exposures 
over multiple epochs. However, the use of type la SNe 
for cosmology i n st rongly lensed regions is fraught with 
difficulty (see p^ . While a failure of ACS would be 
a significant blow to the survey's dark energy science 
objectives, all other components of the CLASH science 
program, including our prime objective of studying the 
distribution of mass profile properties in clusters, would 
be able to be pursued with WFC3 alone. 

If ^CS" fails permanently, we would abandon the dual 
orientation strategy and would, most likely, abandon the 
multiple epoch exposures, allowing the observations of 
each cluster to be completed on a much shorter time- 
frame. By observing each cluster at a single roll angle, 
we would slightly increase the area in which we obtain 
full 16-filter coverage. All of the ACS filters are available 
with WFC3/UVIS, which would allow us to continue our 
16-filter observations with WFC3 alone. However there 
would be noticeable reductions in signal-to-noise in the 
redder fihers (F775W, F814W, F850LP). In order to con- 
tinue to achieve the same depths, we would require ex- 
posure times in the redder filters to be increased by 80%, 
adding about four orbits of integration time to each clus- 
ter. 

The WFC3/UVIS and WFC3/IR fields of view are 
large enough to contain and yield robust photometric 
redshifts for most of the strongly lensed arcs. We would 
lose the extra areal coverage provided by the ACS ob- 
servations, slightly reducing the area of overlap with our 
Subaru images. However, we estimate the impact on our 
ability to measure the mass profiles over the full radial 
range will be negligible. 

5. THE CLASH DATA PIPELINE 

The bulk of our HST data analyses make use of mo- 
saics of globally aligned and co-added images. To accom- 



moer et al.| 


2002 


its input tr 


le cal 



The MosaicDrizzle pipeline takes as 
ibrated FLT files which have been pro- 
duced by calacs and calwfS. The ACS/WFC data, 
however, are first corrected for bias-striping a nd charge 



Koekemoer et al. 2011) supernova programs are tightly 
coordinated and, in tact, share a common pool of re- 
serve orbits from which both programs can draw upon 
for follow-up. This common reserve has been allocated 



transfer efficien cy (CTE) degradation effects (Anderson 
& Bedin 20101. The CTE correction approximately re- 
verses the ettects of charge being trapped and trailed 
across the image during data readout. However, a simi- 
lar procedure has yet to be developed and implemented 
for WFC3/UVIS images. We find that this uncorrected 
CTE can most significantly affect our UVIS photometry 
as follows. Trails from cosmic rays can leak into photo- 
metric apertures of non-detections, artificially boosting 
their observed fiuxes. Objects otherwise expected to be 
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Figure 14. Robustly distinguishing between higli-redshift objects (left) and lower redshift reddened galaxies (right) with the CLASH filter 
set. Photometry in the 16 CLASH filters is plo tted versus w avelength for two galaxies behind CLASH clusters. Left: a spectroscopically 
confirmed z = 6.027 galaxy behind AbcU 383 ( Richard et al.||2011 ) and best fit photometric redshift z = 6.06. Right: a reddened galaxy 
behind Abell 2261 with photo-2 ~ 1.7. Red X's show l-cr detection limits based on the measured flux uncertainties, and blue boxes are 
SED-predicted magnitudes as integrated through each filter. Note the high redshift galaxy SED is roughly fiat in the NIR, while the lower 
redshift SED continues to rise as a function of wavelength through the NIR. This difference, highlighted clearly in the NIR photometry, 
enables us to robustly distinguish high redshift galaxies from lower redshift interlopers. 



UVIS dropouts, or non-detections (based on ACS+IR 
photometry), may commonly have significant detections 
in one or more UVIS filters. We find this can be greatly 
mitigated by adopting a more aggressive rejection of cos- 
mic rays and their trails. The result is that most ex- 
pected UVIS dropouts do actually drop out, though with 
some significant (as high as ~ 5-cr) detections remaining. 

MosaicDrizzle then carries out a sequence of steps 
aimed at aligning the exposures in all the different cam- 
era/filter combinations within each visit and also across 
visits, using a combination of catalog matching and cross- 
correlation. The catalog-matching is done using a deep 
ground-based catalog (typically from Subaru's Suprime 
Cam) as an initial reference point, followed by subse- 
quent matching using catalogs from the CLASH HST 
images themselves, solving for shifts and rotations, which 
yields an accuracy of 0.1-0.2 pixel. Shifts are further re- 
fined using cross-correlation which can achieve accuracies 
to the level of 0.02-0.05 pixel, thus 1-2 mas for ACS im- 
ages, limited essentially at that point by the stability of 
the PSF from one exposure to the next, as well as the ac- 
curacy of the distortion correction which is also removed 
from the images at this step. 

The iterative shift refinement is then used to enable 
cosmic ray rejection to be carried out for all the expo- 
sures in a given filter across multiple epochs, leading to 
a final set of cosmic ray and bad pixel masks. These 
are then weighted according to the sky level in each in- 
put exposure (as modulated by the flatfield / detected 
quantum efficiency variation across the detector for each 
filter), along with the readnoise and accumulated dark 
current, to form an inverse variance image for each ex- 
posure. The drizzle combination is then carried out us- 
ing the inverse variance images as weights, using a square 
kernel as well as a"pixfrac" parameter typically set to 0.8 
which is chosen to be appropriate to the output image 
pixel scale and the number of dithers per filter. An ad- 
vantage of this approach is that it involves only a single 
geometric operation, transforming the pixel values from 
the FLT files directly onto the output frame and avoiding 
any additional convolution. 

It is worth noting that the inverse variance describes 



the expected noise in the absence of any correlated noise; 
in practice the images do have some amount of correlated 
noise (due to the PSF being sampled by the detector pix- 
els to begin with, as well as the single-step transforma- 
tion onto the output image plane), but typically this is 
no more than ~ 10 — 15% of the noise that would be ex- 
pected in the absence of any correlation, and cannot be 
reduced further without making the output image pix- 
els larger. For each cluster, two sets of image mosaics 
are generated — one drizzled to a 30 mas/pixel grid and 
one drizzled to a 65 mas/pixel grid (for both ACS and 
WFC3). We register the images with north up, and the 
images are aligned so that the objects are sampled by 
the same pixels for all the filters. 

We also run a redundant image processing pipeline that 
is built upon th e framework developed for the ACS GTO 
pipehne APSIS ( [Blakeslee et al ||2003[ ) . For the CLASH 
application, APSIS was modified to handle the WFC3 de- 
tectors. The use of this second pipeline provides a useful 
verification of the photometry and astrometry as well 
as a redundant processing facility. The main difference 
between MosaicDrizzle and APSIS is that APSIS does 
not require an external astrometric reference catalog and 
can derive the relative image shifts from the HST data 
alone. To date, both pipelines appear to generate data 
of roughly similar quality with MosaicDrizzle produc- 
ing slightly better alignment precision. 

The SN detection pipeline has been developed jointly 
by the CLASH and CANDELS teams. It wiU be de- 
scribed in more detail in a separate paper. Briefly, the 
SN detection pipeline uses an image-differencing scheme 
to identify SN candidates. Extensive work has already 
been devoted to minimizing false positive detections. IR 
persistence is one potential source for false positives in 
both the SN and high-z galaxy detection projects. To 
address this, a dark exposure is obtained in the Earth- 
occultation just preceding the first visit for all CLASH 
and CANDELS observations. This enables a mask to be 
generated that can flag any suspect pixels in the initial 
exposure of a CLASH visit. Subsequent orbits in the 
visit are flagged based on our actual data. 
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5.1. Object Detection and Characterization 
SExtractor ( [Bertin fc Amouts||1996 1 is used to detect 



or that are only partially detected, we construct manual 
photometric apertures. We force SExtractor to a dopt 



objects and measure their photometry We run SExtrac 
tor (version 2.5.0) in dual image mode, using a detection 
image created from a weighted sum of the ACS/WFC 
and WFC3/IR images. The weights come from the in- 
verse variance images produced by Mosaicdrizzle. We 
do not use the WFC3/UVIS images in the construction 
of the detection image but, of course, run SExtractor 
on the UVIS data to compute source photometry. We 
also create a detection image solely from the WFC3/IR 
images to optimize the search for high redshift (z > 6) 
objects. 

For object detection in the ACS+IR images, we require 
a minimum of 9 contiguous pixels at the level of the ob- 
served background RMS or higher. The detection phase 
background sky level is computed in 5 x 5 grids of cells, 
with each cell being 128 x 128 pixels. For photometry, 
the local sky is estimated from a 24-pixel wide rectangu- 
lar annulus around each detected object. The deblender 
minimum contrast ratio and number of threshold levels 
are 0.0015 and 32, respectively. These parameters were 
chosen after a systematic investigation of possible values 
and yields reasonable performance in minimizing spuri- 
ous detections and suppressing over-deblending of bright 
objects, while achieving reasonable completeness in de- 
tecting faint sources behind bright cluster galaxies. The 
segmentation map for the cluster MACS1149. 6-1-2223 us- 
ing the above parameters is shown in Figure |15[ 

Our image processing effectively rejects cosmic rays 
from the central regions of our UVIS and ACS images 
where we have four or more overlapping exposures. How- 
ever such rejection is not possible in the corners and edges 
of our images, where we have fewer exposures due to our 
observing strategy (two roll angles and dithering). As 
each cosmic ray only affects a single observation, it will 
often only be detected in a single filter. We prune these 
detections from our catalog by rejecting any object with 
only a single B-a detection in one UVIS/ ACS filter, as 
measured by SExtractor. We also reject any object with- 
out any 5-a detections. 

We perform a separate ZR-based detection with more 
aggressive deblending (64 levels of 0.0001 minimum con- 
trast) and background subtraction (3x3 grids of 64 x 64 
pixels). This detection is slightly more sensitive to red- 
der objects, including those at high-redshift. It also per- 
forms slightly better at deblending these fainter objects, 
including those at high-z as well as arcs (strongly-lensed 
galaxies), from brighter nearby cluster galaxies. As all 
WFC3/IR observations are obtained with multiple read- 
outs and cosmic ray rejection, and to preserve our sensi- 
tivity to faint high-redshift objects, we do not prune this 
catalog based on b-a detections. 

Automated detection of strongly lensed galaxies is of- 
ten challenging due to crowding in the cluster core, es- 
pecially for low-surface brightness galaxies stret ched into 



long, thin arcs. P rogress has been made (e.g., Seidel & 
Bartelmann 2007), and we are testing this method tor 



use with CLASH data. Currently, manual intervention 
is still required for those arcs which elude detection. We 
compare the initial source list with arc candidates identi- 
fied both visually and based on our strong lens modeling. 
For arcs that are either missing from the detection list 



these apert ures using the software package SExSeg (Coe 



et al. 2006 1 . Arc photometry is then derived from im- 
ages m which light from the BCG has been modeled and 
subtracted. 

Isophotal apertures are used as they ha ve been shown 
to yield robust co lors ( Bem'tez et al.|2004| ). In the source 
list presented in ^5.3 no aperture corrections have been 
applied to t he magn tudes. The PSF FWHM's span 



20101. The photometric corrections tor 



0.07" -0.15" (|Ford et al.|2003||Sirianni et al.|2003 



Dressel 
variation 

are small for extended objects for all but the faintest 
sources, in part because we use relatively large isophotal 

apertures. 

Ext i nction corrections are derived from the |Schlegel] 
et al.| ( |1998[ ) IR dust emission maps, and the resulting 
A\ coetticients, in AB mag per unit E(B-V), are given 
in Table [5j Flux uncertainties are derived by SExtractor 
using the inverse variance images. We have compared 
the flux errors so derived with the predicted uncertainties 
from the ACS and WFC3 exposure time calculator and 
find excellent agreement. 

5.2. Photometric Redshifts 

We derive photometric re dshift estimates using two in- 
dependent packages: BPZ (lBem'tez"2000^ 'Bem'tez et al 



2004 



Coe et al. 2006) and LePHARE (llbert ct ai. 2006" 



hereafter LPZ). Both software packages use x'^ minimiza- 
tion and template fitting but differ in their specific tem 
plates and their assumed priors. LPZ uses the SED li 



brary optimized for the COSMOS survey (llbert et al. 
2009 1 without template inter polation. BPZ currently uses 



PEU ASE SED templates (Fioc & Rocca-Volmerange 
19971 which have been heavily recalibrated based on 
the FIREWORKS spe ctroscopic and photometric cata- 
log Wuyts et al. ( 2008[ ). BPZ allows for interpolation be- 
tween adjacent templaxes and uses an empirically-derived 
prior on redshift and type based on observed magnitude. 
A full analysis of the accuracy and precision of our pho- 
tometric redshift e stimates will b e discussed in an up- 
coming paper (Jouvel et al.||2012 |. The photoz accuracy 
analyses will be based on over 1, 000 redshifts, obtained 
from the facilities described in ^|6j 

5.3. Source List for MACS1149.6+2223 

We present here an initial version of a represen- 
tative CLASH source list for the z = 0.544 cluster 
MACS1149.6-I-2223. The full hst is available on-line 
(at MAST treasury program archive website) and has 
over 100 parameters per object. In Table W\ we present 
just a small subset of the full list: 100 galaxies within 
1 arcminute of the cluster center. This source list is 



based on the ACS-I-WFC3/IR detection image (see { 5.1 ). 
Column 1 gives the SExtractor object ID, columns 2 
and 3 give the equatorial coordinates (J2000), column 
4 gives the isophotal area of the detection (in 0.065" pix- 
els), column 5 gives the ellipticity of the source as mea- 
sured by SExtractor, column 6 gives the number of pass- 
bands in which the object is detected at a significance 
of 5cr or greater, and column 7 gives the total number 
of bands in which the object is detected. In the case of 
MACS1149.6-f2223, there are 17 HST passbands avail- 
able because the archival imaging included ACS F555W. 
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Figure 15. SExtractor segmentation map for the central region of the cluster MACS1149. 6+2223. Detection of additional objects behind 
the cluster is accomplished by subtracting the best fits to the surface brightness distributions of the brightest of the cluster galaxies. 



Note that even bright sources will sometimes not be de- 
tected in all bands because the footprints of each detec- 
tor on the sky are different (see Figure 111. Columns 
8-10 give the BPZ photo metric redshift estimate, the 



BPZ Odds parameter ( Benftez et al 
measurement presented in Coe et al. 



20041), and the 
12DD6|. The 95% 



confidence limits on the range of the photometric redshift 
estimate is given in the row just below the best fit red- 
shift value. This range is an indicator of the uncertainty 
in the redshift value. In general, BPZ photo-z estimates 
with Odds > 0.90 are the most reliable estimates. Fig- 
ure [16] shows the histogram of the BPZ photometric red- 
shifts for 232 galaxies in the cluster MACSl 149.6-^2223 
with BPZ Odds > 0.90, F850LP magnitude < 26.5 and 
Zphot ^ 3. The histogram shows the cluster peak very 
clearly. 

The remaining columns in Table [7] give the the isopho- 
tal magnitudes (corrected for Galactic extinction) for 
five of the CLASH passbands: F275W, F390W, F850LP, 
F125W, and F160W. The photometric error in the 



isophotal magnitude is given, in parentheses, in the row 
directly below the magnitude value. Photometry has 
been corrected for Galactic extinction using an E(B-V) 
— 0.02297 and all magnitudes given are on the AB pho- 
tometric system (see Table [5] for extinction coefficients 
used) . A magnitude value of 99 indicates a non-detection 
(negative flux) in that band. In this case, the magnitude 
error value gives the la detection limit in that band. A 
magnitude value of —99 indicates the source lies outside 
of the detector field of view or lies within a gap between 
detectors. Source lists for each CLASH cluster are re- 
leased, via MAST, ~ 6 months after the final observation 
for that cluster is acquired. 

6. SUPPORTING OBSERVATIONS 

As discussed above, having both weak and strong lens- 
ing information as well as information about the cluster 
baryonic mass distribution are critical for deriving robust 
mass profiles and concentrations. A partial list of the 
supporting multi-wavelength imaging being used in the 
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Table 7 

Sample of CLASH Source List for MACS1149.6+2223 









Area 








zbpz 






F275W 


F390W 


F850LP 


F125W 


F160W 


ID 


aj20oo 


l5j2000 


(Pixels) 


Ell 


Nsig5 


Nfobs 


[min,max] 


odds 




(err) 


(err) 


(err) 


(err) 


(err) 


4756 


177.4090 


22.4002 


9 


0.36 


4 


17 


2.990 
[2.73,3.09] 


0.63 


1.32 


29.64 
( 1.50) 


29.05 
( 0.65) 


27.18 
( 0.22) 


28.65 
( 0.52) 


28.10 
( 0.30) 


4761 


177.3900 


22.3999 


65 


0.47 


12 


17 


1.701 
[1.62,1.75] 


0.94 


2.44 


26.85 
( 0.49) 


26.73 
( 0.24) 


26.87 
( 0.39) 


25.42 
( 0.08) 


25.65 
( 0.08) 


4762 


177.3900 


ZZ.4UUZ 


10 


0.45 


4 


17 


1.717 
[1.49,1.77] 


0.59 


0.63 


29.20 
( 1.53) 


28.16 
( 0.48) 


27.16 
( 0.33) 


26.82 
( 0.17) 


27.11 
( 0.19) 


4763 


177.3900 


on /i nnn 

22.4000 


8 


0.42 


1 


17 


0.765 
[0.76,2.22] 


0.12 


0.63 


99.00 
(28.12) 


28.19 
( 0.44) 


27.28 
( 0.32) 


27.26 
( 0.21) 


27.36 
( 0.20) 


4764 


177.3910 


22.4002 


29 


0.27 


12 


17 


1.222 
[1.12,1.32] 


0.63 


1.98 


99.00 
(28.01) 


27.69 
( 0.34) 


26.44 
( 0.18) 


26.24 
( 0.10) 


26.26 
( 0.09) 


4765 


177.4080 


nn ^ nn o 

22.4002 


15 


0.16 


8 


17 


0.745 
[0.67,4.03] 


0.31 


1.47 


99.00 
(28.04) 


99.00 
(28.89) 


26.59 
( 0.19) 


26.94 
( 0.18) 


26.85 
( 0.14) 


4773 


177.4140 


22.4002 


4 


0.31 


1 


17 


1.150 
[1.25,2.41] 


0.14 


0.87 


28.89 
( 0.95) 


27.87 
( 0.25) 


27.97 
( 0.39) 


28.95 
( 0.60) 


28.03 
( 0.27) 


4778 


177.3880 


22.3999 


152 


0.13 


12 


17 


0.453 

[0.43,0.47] 


0.98 


4.00 


25.34 
( 0.20) 


26.62 

( 0.30) 


23.78 
( 0.04) 


23.48 
( 0.02) 


23.30 

( 0.01) 


4780 


177.4130 


oo j4 nnn 

22.4000 


260 


0.03 


13 


17 


0.685 
[0.66,0.72] 


0.91 


5.22 


26.26 
( 0.53) 


26.32 
( 0.29) 


23.35 
( 0.03) 


23.23 
( 0.02) 


23.00 
( 0.01) 


4782 


1 T'7 Q or\c\ 

177.3890 


z2.4Uul 


8 


0.06 


2 


17 


2.344 
[2.20,2.76] 


0.32 


0.28 


99.00 
(28.58) 


28.35 
( 0.36) 


28.67 
( 0.67) 


27.85 
( 0.25) 


27.07 
( 0.11) 


4784 


177.4100 


or\ >i nnn 

22.4000 


10 


0.51 


4 


17 


0.887 
[0.51,1.00] 


0.37 


1.86 


26.50 
( 0.28) 


27.07 
( 0.22) 


26.65 
( 0.22) 


27.06 
( 0.23) 


27.52 
( 0.29) 


4786 


177.3870 


22.4001 


5 


0.03 


2 


17 


3.898 
[3.65,4.08] 


0.64 


1.30 


99.00 
(28.87) 


30.42 
( 1-21) 


28.65 
( 0.53) 


29.30 
( 0.60) 


29.74 
( 0.72) 


4787 


177.4100 


oo onno 

22.3998 


164 


0.53 


12 


17 


0.507 
[0.46,0.54] 


0.80 


3.54 


25.89 
( 0.34) 


27.03 
( 0.43) 


24.17 
( 0.06) 


24.19 
( 0.04) 


24.23 
( 0.04) 


4788 


177.414U 


z2.4Uul 


7 


0.20 


2 


17 


2.472 
[1.93,2.68] 


0.40 


0.83 


99.00 
(28.45) 


27.85 
( 0.27) 


27.84 
( 0.38) 


28.24 
( 0.38) 


28.82 
( 0.51) 


4792 


177.3820 


22.4000 


5 


0.45 


6 


17 


0.587 
[0.45,2.98] 


0.31 


0.67 


99.00 
(27.95) 


28.13 
( 0.49) 


27.33 
( 0.40) 


26.98 
( 0.20) 


26.83 
( 0.15) 


4793 


177.3910 


nn onno 

22.3998 


387 


0.29 


14 


17 


3.248 
[3.23,3.27] 


1.00 


1.93 


25.38 
( 0.28) 


25.03 
( 0.11) 


23.05 
( 0.03) 


23.01 
( 0.02) 


22.69 
( 0.01) 


4794 


177.3910 


on onnn 

22.3999 


26 


0.36 


11 


17 


3.274 
[3.25,3.42] 


0.92 


1.34 


99.00 
(28.36) 


28.92 
( 0.64) 


26.96 
( 0.20) 


26.71 
( 0.12) 


26.44 
( 0.08) 


4799 


177.3830 


nn onno 

22.3998 


115 


0.21 


16 


17 


1.140 
[1.13,1.17] 


1.00 


2.94 


25.07 
( 0.15) 


24.97 
( 0.07) 


24.06 
( 0.05) 


23.98 
( 0.03) 


23.74 
( 0.02) 


4800 


177.3840 


nn onno 

22.3998 


132 


0.20 


13 


17 


0.714 
[0.67,0.74] 


0.88 


2.33 


25.64 
( 0.22) 


26.46 
( 0.23) 


24.59 
( 0.07) 


24.71 
( 0.05) 


24.78 
( 0.05) 


4801 


177.3840 


nn onn'^ 

22.3997 


472 


0.10 


15 


17 


0.528 
[0.46,0.54] 


0.82 


4.85 


24.30 
( 0-14) 


24.56 
( 0.09) 


22.75 
( 0.03) 


22.72 
( 0.02) 


22.63 
( 0.01) 


4805 


177.3950 


nn ,1 nnn 
22.4UUU 


29 


0.21 


12 


17 


0.906 
[0.88,0.94] 


0.95 


3.55 


27.30 
( 0.49) 


27.30 
( 0.26) 


25.84 
( 0-12) 


26.34 
( 0.12) 


26.50 
( 0-12) 


4806 


177.4070 


f^r* ^ nnn 

22.4000 


58 


0.34 


14 


17 


0.776 
[0.74,0.84] 


0.81 


2.52 


26.09 
( 0.25) 


26.48 
( 0.18) 


25.48 
{ 0-11) 


24.99 
( 0.05) 


24.91 
( 0.04) 


4812 


177.4150 


on onnn 

22.3999 


69 


0.19 


10 


17 


0.347 

[0.30,0.39] 


0.74 


1.53 


99.00 

(27.49) 


27.99 

( 0.61) 


25.82 
( 0.16) 


24.93 

( 0.05) 


24.81 

( 0.04) 


4813 


177.4050 


on onno 

22.3998 


321 


0.04 


14 


17 


0.505 
[0.49,0.52] 


1.00 


2.09 


27.06 
( 0.94) 


25.72 
( 0.19) 


22.59 
( 0.02) 


22.17 
( 0.01) 


21.85 
( 0.01) 


4820 


177.3840 


z2.4UU7 


76 


0.02 


11 


17 


0.772 
[0.71,0.84] 


0.73 


3.08 


99.00 

(27.36) 


26.97 
( 0.32) 


25.31 
( 0.12) 


25.03 
( 0.06) 


25.00 
( 0.05) 


4821 


177.3920 


no onnr' 

22.3995 


76 


0.11 


10 


17 


0.742 
[0.54,0.81] 


0.49 


2.22 


27.21 
( 0.68) 


99.00 
(28.17) 


25.56 
( 0.14) 


24.99 
( 0.06) 


24.78 
( 0.04) 


4822 


177.3920 


on onno 

22.3998 


32 


0.34 


9 


17 


0.560 
[0.34,3.09] 


0.11 


0.89 


27.40 
( 0.69) 


27.48 
{ 0.41) 


26.13 
( 0.21) 


26.52 
( 0.19) 


26.01 
( 0.11) 


4827 


177.3860 


on onno 

22.3998 


5 


0.09 


1 


17 


0.650 
[0.53,0.82] 


0.40 


0.71 


27.72 
( 0.37) 


28.66 
( 0.42) 


28.34 
( 0.49) 


28.44 
( 0.36) 


28.33 
( 0.29) 


4832 


177.4120 


oo onn£? 

22.3996 


66 


0.54 


14 


17 


2.739 
[2.64,2.75] 


0.96 


1.41 


99.00 
(27.56) 


26.10 
( 0.13) 


25.33 
( 0.10) 


25.31 
( 0.07) 


25.03 
( 0.05) 


4841 


177.4040 


22.3996 


74 


0.37 


14 


17 


0.371 
[0.34,0.47] 


0.70 


5.95 


26.80 
( 0.57) 


27.16 
( 0-41) 


24.63 
( 0.07) 


25.42 
( 0.10) 


25.07 
( 0.06) 


4843 


177.3890 


22.3996 


21 


0.55 


6 


17 


0.768 
[0.64,1.03] 


0.37 


0.92 


99.00 
(27.85) 


27.84 
( 0.42) 


26.57 
( 0.23) 


26.19 
( 0.11) 


26.03 
( 0.08) 


4853 


177.4140 


22.3994 


56 


0.28 


13 


17 


2.201 
[2.19,2.33] 


0.82 


2.38 


27.57 
( 0.75) 


25.83 
( 0.10) 


26.66 
( 0.29) 


25.94 
( 0.11) 


25.27 
( 0.05) 


4856 


177.4050 


22.4009 


21 


0.33 


8 


17 


0.767 
[0.63,0.92] 


0.41 


2.30 


26.01 

( 0.20) 


27.68 
( 0.40) 


26.15 
( 0.17) 


26.34 
( 0.14) 


26.01 
( 0.09) 


4868 


177.4060 


22.3994 


21 


0.17 


7 


17 


0.558 
[0.51,0.97] 


0.27 


1.64 


27.33 
( 0.53) 


27.98 
( 0.47) 


26.42 
( 0.19) 


26.68 
( 0.17) 


26.26 
( 0.10) 


4869 


177.4060 


22.3991 


215 


0.03 


12 


17 


0.455 
[0.43,0.47] 


0.99 


2.13 


25.70 
( 0.33) 


26.90 
( 0.44) 


23.48 
( 0.04) 


23.16 
( 0.02) 


22.90 
( 0.01) 



Note. — The full source list with over 100 parameters per source is available on-line from lOP and from the MAST CLASH Treasury Archive 

webpage at http://archive.stsci.edu/prepds/clash/ 



All magnitudes are on AB system. 
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1 2 3 

BPZ Photometric Redshift 

Figure 16. Histogram of tlie BPZ pliotometric redsliifts 
for galaxies witliin 1 arcminute of tlie center of the cluster 
MACS1149. 6+2223. The 232 galaxies included here have BPZ 
Odds > 0.90, F850LP magnitude < 26.5 and z„hot < 3. The 
histogram shows the cluster peak very clearly. The mean spec- 
troscopic redshift of the cluster (z = 0.544) is denoted by the 
vertical red line. The average photometric redshift in the range 
0.4 < Zphot < 0.65 is {zph) = 0.515. The accuracy of our photo- 
metric redshifts is exp ected to improve fur ther as our spectroscopic 
sample increases. See jjouvel et alT| | |2012[ l for details. 

CLASH program is summarized in Table |8| All CLASH 
clusters have X-ray imaging as well as wide-field multi- 
band ground-based opti cal imaging. The X-ray imaging 
is from Chandra/ ACIS (Garmire et al. 2003) and some 
of the clusters have XMM/EPIC pruder eraL]|200T 
Turner et aL][200l| imaging as well 



I'he ground- based optical imaging is, for 24 out of 
the 25 clusters, from the Subaru Observatory (and, in 
particular, performed with the SuprimeCam instrument 
( Miyazaki et al.|1998[ )). The Subaru SuprimeCam imag- 
ing is currently available in three or more filters for 22 
of the 24 clusters reachable from Mauna Kea and we are 
pursuing additional optical filter coverage for the two re- 
maining clusters. Wide-field optical imaging for AbcU 
2261 from Subaru, for instance, was supplemented with 
SDSS iz band imaging from the Kitt Peak 4-m Mayall 
telescope using the Mosaicl camera. RXJ2248. 7-4431, is 
too far south to reach with Subaru but has multi-band 
imaging from the wide-field imager (WFI) on the Euro- 
pean Southern Observatory's 2.2m telescope. The multi- 
band optical imaging enables robust selectio n of back- 
ground galaxies for weak lensing analysis (e.g., Medezin- 
ski et al. 2007). 

We have also acquired mm-wave imaging to map the 
Sunyaev-Zel'dovich Effect (SZE) in order to measure the 
integrated line-of-sight gas pressure towards each cluster. 
The SZE data come primarily from two facilities - the 
Boloca m instrument at the Caltech Submillime ter Obser- 
vatory (iGolwala et al.|2009[ |Sayers et al.|2011j ) and Mul- 
tiplexec TSQUlD/ I'ES T^rray at INinety Gigahertz (MUS- 
TANG ; [Dicker et al.| [2008, 2009) at the Robert C. Byrd 
Green Bank Telescope. All 25 CLASH clusters have a de- 
tected SZE signal with the lowest SNR^ 7 and most with 
SNR > 10. The combination of the X-ray and mm-wave 
observations allow the mass scaling relations to be ac- 



curately calibrated f or use in cosmological surveys (e.g., 
Okabe et aL]|2010b[ ). The X-ray and SZE data will be 
used to measure gas density profiles for subtraction from 
lensing-deriv ed total mass pro files to yield DM-only mass 
profiles (e.g., Lemze et al. 2008). These data also allow 
assumptions ot hydrostatic equilibrium to be tested and 
constraints on thermal motion to be made from measure- 



2010 



Molnar et al. 



Bpitzer Space ite. 



2010) 



ments of excess pressure ( Umetsu et al. 2009 Zhang et al. 



escope (Werner et al. 2004) data is 
available for most of the CLAsH clusters trom the IRAC 
Lensing Survey (E. Egami, P.I.). The few remaining clus- 
ters that lacked sufhciently deep Spitzer data will be ob- 
served in our cycle 8 ICLASH program (R. Bouwens, 
P.I.). The Spitzer data are critical for characterizing 



the hi gh-z galaxy pop ulation (e.g., see Richard et al 



pOTTi). Deep IRAC ( [Fazio et alT]|2004[ ) imaging sam 
pies tlie rest-frame optical in high- z galaxies, allowing us 
to perform stellar population modeling and derive stel- 
lar masses, and enables robust discrimination between 
star- forming z > 7 galaxies and dusty/evolved galaxies 
at z ~ 2. 

We have also been awarded 225 hours of time on the 
VLT as part of a large program (P. Rosati, P.I.) to ob- 
tain VIMOS spectroscopy of the 14 southern clusters 
(^ < -1-1.5°). In addition, we are acquiring cluster and 
lensed galaxy redshifts from spectrographs on Magellan 
(GISMO; LDSS), LBT (MODS), MMT (Hectospec), and 
Palomar (double spectrograph). 

7. SUMMARY 

The Cluster Lensing And Supernova survey with Hub- 
ble will produce a major advance in our understand- 
ing of the DM power spectrum and the internal struc- 
ture of cluster halos on scales from 10 kpc to 2 Mpc. 
The precision to which these measurement are being 
made will provide an unprecedented foil against which we 
will challenge and ultimately expand our current ideas 
about structure formation and the nature of dark en- 
ergy. In November 2010, the CLASH MCT program 
initiated its three year observing plan to obtain deep 
(20-orbit) 16-band HST imaging for each of 25 galaxy 
cluster cores. The clusters in our sample are massive 
(5 X 10" < M„r/MQ < 3 X 10^5), span a range of (in- 
termediate) redshifts (0.18 < z < 0.9), and most (20 
of the 25) were selected based on their X-ray properties 
(reasonably relaxed and > 5 keV). The 16-band HST 
imaging yields precise (2%(1 -I- z)) photometric redshifts 
for all galaxies brighter that F775W AB mag 26, includ- 
ing hundreds of strongly lensed galaxies. When com- 
bined with weak lensing maps from wide field Subaru 
imaging and constraints on the baryonic gas distribution 
from mm-wave and X-ray imaging, CLASH data will al- 
low the cluster mass profiles to be tightly constrained 
over a large range of scales, yielding robust measure- 
ments of their central density concentrations. We also 
expect to provide important constraints on the degree of 
substructure in the DM distribution by studying galaxy- 
scale lensing. The strong lens magnifying power of our 
clusters should also enable detection of dozens of rela- 
tively bright (m < 26 AB) z > 7 galaxies, including some 
bright enough for spectroscopic follow-up with ground- 
based telescopes and, certainly, many that can be studied 
with the spectrographs on JWST. Parallel observations 
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Table 8 

Partial List of Supporting Observations for CLASH at Time of Publication 



















Suprime Cam 


Suprime Cam 


Chandra 


XMM 


Spitzer^ 


Bolocani 


Cluster 


(Mpc) 


Filters 


(ksec) 


(ksec) 


(ksec) 


(ksec) 


X-ray Selected Clusters: 














Abell 209 


5.68 


BVRIZ 


20 




(31.2)'' 


66 


Abell 383 


5.26 


BVRIZ 


50 




20.4 


91 


MACS0329.7-0211 


9.68 


BVRZ 


70 




(31.2)'' 


39 


MACS0429. 6-0253 


9.30 


VRI 


24 




(31.2)'' 


62 


MACS0744.9+3927 


11.90 


BVRIZ 


90 


138 


18.0 


62 


Abell 611 


7.28 


BVRIZ 


40 


28 


18.0 


79 


MACS1115. 9+0129 


8.34 


BVRIZ 


56 




(31.2)'' 


56 


Abell 1423 


5.82 


VI 


36 


31 




TBA'' 




y.uu 


BVRIZ 


24 






41 


CLJ1226.9-(-3332 


13.04 


BVZ 


75 


120 


17.0 


68 


MACS1311.0-0310 


10.20 


R 


85 




(31.2)'' 


41 


rtAJ iv>4/ .0-ii4o 


9.68 




185 


36 


19.2 


59 


iVi jt. o 1 4 z o . o -r 4 u 4 


iU. ( z 


Jj V IXIZj 


1 /I n 

i'-tU 






oo 






Jj V ixIZj 


A n 




on o 


Oo 


A/T A r^QI ion Q 1 QC^Qfi 
iVi/AL^oi r ZU.o-hoOoD 


Q on 
o.yu 


Jj V xxIZj 


fin 
DU 




on n 


Do 


Abell zzDl 


O.UO 


r> V rt 


oo 


QQ 
OO 


oo /I 


4D 


MACsl931. 8-2635 


8.34 


BVRIZ 


114 




(31.2)'' 


63 


RXJ2129. 7-1-0005 


6.26 


BVRIZ 


40 




18.2 


49 


MS2137-2353 


7.70 


BVRIZ 


150 




20.4 


47 


RXJ2248. 7-4431 (Abell 1063S) 




Not accessible 


26 




20.4 


20 






See note c 










High Magnification Clusters: 














MACS0416.1-2403 


9.30 


BRZ 


16 




(31.2)'' 


29 


MACS0647.8+7015 


11.10 


BVRIZ 


40 




18.0 


43 


MACS0717.5+3745 


10.76 


BVRIZ 


104 




21.4 


47 


MACS1149.6+2223 


10.72 


BVRIZ 


20 




18.0 


45 


MACS2129.4-0741 


10.96 


BVRIZ 


40 




19.5 


47 



Note. — Observations with MUSTANG (NRAO/GBT) pending. Spectroscopic observations summarized in text. 
^ Times shown are the total for both IRAC channels 1,2 (3.6 and 4.5 microns). 
^ These 8 clusters will be observed with the indicated exposures in Spitzer cycle 8. 

RXJ2248. 7-4431 has UBVRIZ imaging with the WFI on the 2.2-meter telescope at La Silla. Total exposure times 
in each band are 10-11 hours. 

SZE data for A1423 will bo acquired in spring 2012. 



may detect up to ~ 30 z > 1 Type la supernovae that, 
when combined with z > 1 SNe la detected by the CAN- 
DELS program, will provide new constraints on the time- 
dependence of w and a dramatic improvement in our un- 
derstanding of SNe la evolution in the matter-dominated 
universe (z > 1.5). 

CLASH data will also provide the mass-calibrators for 
the next generation of big cosmological surveys such as 
the Dark Energy Survey (DES), Sunayev-Zeldovich sur- 
veys (e.g., South Pole Telescope) and next generation X- 
ray cluster surveys. The DES will measure the evolution 
of the space density of clusters (among other probes) as 
a way to measure Dark Energy. These cluster counts tell 
us much about cosmological parameters through their 
impact on both the volume and the growth of perturba- 
tions. To do so, however, requires that the cluster mass 
scaling relations be well calibrated. The CLASH mea- 
surements of the enclosed mass in each cluster will be 
invaluable in calibrating the relation between the cluster 
mass and a number of observable mass proxies. 

MCT programs represent a very large investment of 
HST observing time. These data are intended to be a 
community resource, and as such they have no propri- 
etary period. The zero proprietary period policy also 
applies to our cycle 8 Spitzer observations. All HST and 
Spitzer images may be downloaded immediately after 



they are obtained via the MAST and IRSA archives, re- 
spectively. The CLASH collaboration is producing high- 
level science products that will be publicl y d istributed 
via mast's Treasury Program archive sitep^ This pa- 
per includes the first public release of a CLASH source 
fist for the cluster MACS1149.6+2223 as well as uni- 
formly re-derived X-ray properties for the 25 clusters in 
the sample. The CLASH HLSP include co-added and 
mosaicked images, weight maps, inverse variance maps, 
high-resolution color images, and source lists. At the 
end of the survey we will also release the correspond- 
ing registered X-ray surface brightness images, co-added 
Spitzer Space Telescope images, co-added SZE data, the 
calibrated lensing-derived mass and source magnification 
maps, the source lists with all photometric redshift infor- 
mation, SN coordinates and their light curves and grism 
spectra, and the associated spectroscopic redshift cata- 
logs. The CLASH data provide a vast legacy archive for 
studies of the formation and evolution of cosmic struc- 
ture. 

We are especially grateful to our program coordinator 
Beth Perrillo for her expert assistance in implementing 
the HST observations in this program. We thank Jay 

^® http:/ /archive. stsci.edu/hst/tall. html 
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Anderson and Norman Grogin for providing the ACS 
CTE and bias striping correction algorithms used in our 
data pipeline. Finally, we are indebted to the hundreds 
of people who have labored many years to plan, develop, 
manufacture, install, repair, and calibrate the WFC3 and 
ACS instruments as well as to all those who maintain and 
operate the Hubble Space Telescope. 

The CLASH Muhi-Cycle Treasury Program (GO- 
12065) is based on observations made with the 
NASA/ESA Hubble Space Telescope. The Space Tele- 
scope Science Institute is operated by the Association 
of Universities for Research in Astronomy, Inc. under 
NASA contract NAS 5-26555. ACS was developed un- 
der NASA contract NAS 5-32864. 
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AYA2010-22111-C03-00, funding from the Junta de An- 
daluca Proyecto de Excelencia NBL2003, INAF con- 
tracts ASI-INAF 1/009/10/0, ASLINAF 1/023/05/0, 
ASI-INAF 1/088/06/0, PRIN INAF 2009, and PRIN 
INAF 2010, NSF CAREER grant AST-0847157, the 
UK's STFC, the Royal Society, the Wolfson Foundation, 
and National Science Council of Taiwan grant NSC97- 
2112-M-001-020-MY3. AZ acknowledges support by the 
John Bahcall excellence prize. LI acknowledges support 
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Origin and Structure of the Universe program. 
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